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Chapter  1 

Introduction  and  organization 


This  Final  Technical  Report  summarizes  the  main  results  of  the  work  performed  at  Rutgers 
under  Grant  N00014-94-1-0854. 

1.1  Organization  of  the  Report 

The  main  technical  accomphshments  are  outhned  in  in  Chapter  2,  first  in  summary  form,  then 
in  greater  detail. 

All  of  the  work  reported  here  has  been  written  up  in  the  form  of  papers  for  submission  to 
refereed  journeJs.  The  work  reported  in  Chapter  3  below  will  form  the  basis  for  a  paper  to  be 
submitted  for  pubhcation.  Chapters  4  and  5,  have  been  accepted  for  pubhcation  as,  respectively, 

•  M.  A.  Grinfeld  cind  A.  N.  Norris,  “Hamiltonian  and  Onsageristic  approaches  in  the  nonUn- 
ear  theory  of  fluid  permeable  elastic  continua”,  Int  J.  Engng,  Sc, 

•  M.  A.  Grinfeld  and  A.  N.  Norris,  “Acoustoelasticity  theory  and  apphcations  for  fluid- 
saturated  porous  media”,  J,  Acoust,  Soc,  Am, 

Chapter  6  has  been  submitted  for  publication,  as 

•  M.  A.  Grinfeld  and  A.  N.  Norris,  “Waves  ctnd  flutter  in  fluid  conveying  elastic  media”,  J, 
Fluids  Struct, 

One  paper  has  appeeired  in  print: 

•  A.  N.  Norris  and  M.  A.  Grinfeld,  “Nonlinear  poroelasticity  for  a  layered  medium”,  J, 
Acoust,  Soc,  Am,  98,  1138-1146,  1995. 

A  copy  of  this  is  attached  to  the  report.  Two  other  papers  have  been  completed  on  the  topic  of 
nonlinear  elasticity  of  granular  media,  in  collaboration  with  David  Johnson  of  Schlumberger-DoU 
Research.  They  are  not  included  here,  but  are  cited  in  Section  1.2.2  below. 
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1.2  Personnel,  publications  and  presentations 

1.2.1  Personnel 

1.  PI:  Andrew  N.  Norris 

2.  Dr.  Michael  Grinfeld 

Visiting  Research  Professor,  Department  of  Mechanical  &  Aerospace  Engineering 
Rutgers  University 


1.2.2  Publications  resulting  from  Grant  N00014-94-1-0854 

1.  Nonlinear  poroelasticity  for  a  layered  medium,  A.  N.  Norris  and  M.  A.  Grinfeld, 
J.  Acoust.  Soc.  Am.  98,  1138-1146, 1995. 

2.  Acoustoelasticity  theory  and  applications  for  fluid-saturated  porous  me¬ 
dia,  M.  A.  Grinfeld  and  A.  N.  Norris,  J.  Acoust.  Soc.  Am.  {accepted  for  publication). 

3.  Hamiltonian  and  Onsageristic  approaches  in  the  nonlinear  theory  of  fluid 
permeable  elastic  continua,  M.  a.  Grinfeld  and  A.  N.  Norris,  Int.  J.  Engng.  Sci. 
{accepted  for  publication). 

4.  Acoustics  of  fluid  conveying  elastic  structures,  M.  A.  Grinfeld  and  A.  N.  Norris, 
J.  Fluids  Struct,  (submitted). 

5.  Nonlinear  elasticity  of  granular  media,  A.  N.  Norris  and  D.  L.  Johnson,  J.  Appl. 
Mech.  {accepted  for  publication).  [Supported  in  part  by  N00014-94-1-0854.] 

6.  Rough  elastic  spheres  in  contact,  D.  L.  Johnson  and  A.  N.  Norris,  J.  Mech.  Phys. 
Solids  {submitted).  [Supported  in  part  by  N00014-94-1-0854.] 

1.2.3  Presentations  at  conferences 

1.  A.N.  Norris,  Variational  approaches  in  the  nonlinear  theory  of  fluid-permeable 
ELASTIC  CONTINUA,  Acoustical  Society  of  America,  St.  Louis,  MO,  Nov.  1995. 

2.  M.A.  Grinfeld,  Acoustics  of  fluid  conveying  elastic  structures.  Acoustical  Soci¬ 
ety  of  America,  Washington,  D.C.,  June  1995 

3.  M.A.  Grinfeld,  Nonlinear  theory  of  poroelasticity  and  applications.  Third  SIAM 
Conference  on  Mathematical  and  Computational  Issues  in  the  geosciences,  San  Antonio, 
February  1995. 


Chapter  2 

Summary  of  Research 
Accomplishments 


2.1  Abstract  Summary 

Considerable  progress  has  been  made  in  establishing  the  appropriate  master  system  of  equations 
suitable  for  modeling  shock  waves  and  nonlinear  wave  phenomena  in  fluid-permeable  porous 
sohds.  The  major  accomplishments  are  as  follows,  (i)  Numerical  experiments  establishing  the 
existence  of  shock  waves  in  a  ID  system  with  both  Darcy  and  Navier-Stokes  viscosities,  (ii) 
The  general  form  of  the  equations  of  motion  have  been  deduced  from  Hamilton’s  principle  of 
Least  Action  combined  with  Onsager’s  method  of  irreversible  thermodynamics,  (iii)  We  ana¬ 
lyzed  the  specific  case  of  a  layered  fluid/solid  medium  in  depth,  and  obtained  explicit  formulas 
for  the  nonlinejir  straiin  energy  and  the  permeability  operator,  (iv)  We  have  extended  the  meth¬ 
ods  of  acoustoelasticity  to  poroelastic  media,  (v)  We  have  examined  the  stability  of  acoustic 
disturbances  in  a  simple  model  of  a  porous  medium  with  a  mean  flow  through  the  pores. 


2.2  Summary  of  Progress 

The  two  major  objectives  in  our  research  on  shock  waves  in  poroelastic  fluid-filled  media  are: 

1.  To  establish  Ein  appropriate  master  system  (MS)  of  nonlinear  PDEs,  «ind 

2.  Analysis  of  physical  processes  within  the  shock  tremsitional  layer,  and  development  of 
jmalyticcil  and  numerical  tools  for  solving  the  MS. 

A  MS  suitable  for  the  study  of  nonlinear  phenomena  in  fluid  permeable  elastic  media  (2ind  shock 
wave  propagation,  in  particular)  should  satisfy  the  following  principles: 

(a)  It  should  produce  sufficiently  smooth  and  physically  acceptable  finite  solutions  of  initial- 
boundeiry-vcilue  problems  for  eill  t  >  0; 

(b)  Its  derivation  should  be  conceptually  simple  and  rely  on  a  minimal  amount  of  ad  hoc  as¬ 
sumptions; 
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(c)  It  should  describe  all  known  effects  of  fluid-permeable  elastic  media; 

(d)  It  should  include  contemporary  nonlinear  theories  one-component  of  fluids  and  solids  as  par¬ 
ticular  cases. 

Despite  numerous  attempts  of  our  predecessors  in  the  literature  we  were  not  able  to  find  a  MS 
satisfying  requirements  (a)-(d). 

We  believe  that  we  have  made  considerable  progress  in  establishing  the  appropriate  MS  satisfying 
the“a-d”  requirements.  We  now  have  a  reliable  framework  for  studying  a  wide  variety  of  nonlinear 
phenomena  of  wave  propagation  in  fluid-permeable  elastic  media.  The  principal  results  of  our 
work  Eire 

•  Numerical  experiments  establishing  the  existence  of  shock  waves  in  a  ID  system  with  Darcy 
and  Navier-Stokes  viscosities  -  see  Section  2.3.1, 

•  The  general  form  of  the  MS  deduced  from  Hamilton’s  principle  of  Least  Action  combined 
with  Onsager’s  method  of  irreversible  thermodynamics  -  see  Section  2.3.2, 

•  A  specific  form  of  the  MS  for  layered  media  obtained  by  the  method  of  homogenization  - 
2.3.3, 

•  Application  and  extension  of  the  methods  of  acoustoelasticity  to  fluid-filled  poro-elastic 
media  -  see  Section  2.3.4, 

•  A  study  of  the  influence  of  mean  flow  through  a  porous  medium  on  the  stability  of  acoustical 
disturbances.  The  anedysis  shows  that  flutter- type  instabilities  are  possible  -  see  Section 
2.3.5. 

The  final  four  items  have  been  written  up  and  either  submitted  or  published  [1,  2,  3,  4].  The 
first  is  stiU  under  active  investigation,  but  will  be  submitted  in  the  near  future. 


2.3  Detailed  Summary  of  Technical  Accomplishments 

2.3.1  Shock  wave  structure 

The  key  issue  regarding  shock  waves  is  the  question  of  the  shock  wave  structure,  with  the  following 
objectives:  (i)  To  understtind  the  interaction  of  physiceil  processes  within  the  narrow  tramsitional 
layer  of  the  shock,  and  (ii)  to  simplify  the  mathematical  and  numericeil  analysis  by  establishing 
the  analogy  of  Hugoniot  conditions  at  the  shock.  Mathematically,  the  “problem  of  the  shock  wave 
structure”  is  equivalent  to  the  study  of  traveling  waves  with  stationary  profiles.  Because  the  MS’s 
are  rather  complicated  we  first  considered  the  model  equation  (ME)  Ut  +  UUx  +  AU  -k nUxx  =  0. 
The  ME  has  mziny  features  of  the  MS:  it  is  nonlinear  and  possesses  both  Darcy-like  (AU)  and 
Navier-like  (fiUxx)  viscous  effects. 

Our  analysis  of  traveling  wave  solutions  for  the  ME  indicates  (surprisingly!)  that  the  ME  does 
not  permit  smooth  traveling  waves  of  finite  Eimplitude  because  of  the  influence  of  the  Darcy 
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dissipation.  This  finding  rciises  several  conceptual  problems.  Are  there  shock-like  solutions  in 
fluid-permeable  media?  If  “yes”  what  is  the  analogy  of  the  Hugoniot  conditions?  How  shotild  we 
scale  the  parameters  of  the  ME  and  MS  in  order  to  obtain  finite  solutions  of  the  shock- structure- 
problem?  We  decided  to  address  these  problems  with  the  help  of  numericzd  experiments  for  the 
ME.  Our  simulations  clearly  show  the  existence  of  the  shock-like  solutions.  We  are  currently  con¬ 
tinuing  the  numerical  modeling  of  the  model  ID  equation  in  order  to  “experimentally”  establish 
the  Hugoniot  shock  conditions.  For  full  detciils  see  Chapter  3  below. 


2.3.2  Principle  of  Least  Action  for  multi- phase  media 

It  is  common  practice  to  formulate  nonlineeir  equations  of  solids  in  reference  or  Lagrangian  coor¬ 
dinates.  Fluid  dyncimics,  on  the  other  hand,  is  more  naturadly  considered  in  current  or  Eulerian 
coordinates.  The  reasons  are  clear:  in  solids  the  material  particles  adhere  to  one  another,  so 
that  the  reference  description  always  provides  a  continuous  mapping  to  the  current  configuration 
of  peirticles.  Not  so  for  fluids,  where  leirge  scale  shearing  can  separate  neighboring  particles  in 
the  reference  system.  The  same  is  true  of  fluid-solid  composite  media,  because  the  fluid  can 
slide  past  solid  particles.  The  use  of  Eulerian  coordinates  for  such  media  therefore  appears  to 
be  the  most  convenient  procedure  and  maybe  the  only  feasible  one.  Despite  this  obvious  need 
for  2in  EulericOi  description,  the  many  previous  works  on  the  development  of  poroelastic  theories 
Me  very  weak  on  this  point.  They  invariably  obteiin  the  equilibrium  equations  by  postulating 
stresses,  momenta,  and  various  interaction  effects  such  as  a  “momentum  supply”  for  each  phase. 

We  have  taken  a  more  fundaimented  point  of  view:  that  the  equations  governing  the  system 
should  all  follow  unambiguously  from  the  principle  of  least  action  or  Haimilton’s  principle.  Our 
starting  point  is  a  Lagramgian  density  per  unit  mass  of  the  entire  system,  L,  which  is  assumed 
to  be  a  function  of  the  following  quantities  for  each  phase  a  (where  a  can  signify  fluid,  solid,  or 
any  one  of  a  number  of  constituents  in  a  multiphase  medium):  The  density  pa,  the  velocity  Va, 
and  the  displacement  gradient  Vua-  Conservation  of  mass  is  considered  fundamentcd  for  each 
species,  but  no  other  balance  laws  are  invoked  in  our  treatment.  The  principle  of  least  action  is 
based  on  the  stationarity  of  the  “Action”  -  /dt  JdV  L,  and  offers  the  most  general  framework  for 
setting  up  equilibrium  equations  and  boundeuy  conditions.  There  Me  some  formidable  techniced 
difficulties  involved  in  applying  this  principle  in  Euleri^m  vMiables  and  to  multi-phase  media. 
We  have  overcome  these  and  have  found  that  the  general  form  of  the  equilibrium  equations  for 
a  two  component  medium  must  be  of  the  form 

Here,  p  =  pi  +  P2  is  the  tot2il  density,  eind  Pa  is  the  elastic  stress  in  a  which  follows  from 
the  LagrangiJin  density  as  Pjia  =  ^jiPPa^  -  p{hi  -  vector  Ga  is  a  more 

complicated  but  explicit  function  of  the  Lagrangian,  but  it  cein  be  interpreted  as  em  interaction 
force  between  the  two  phases  because  of  the  universed  identity  Gi  -t-  G2  =  0. 
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Equation  (2.1)  is  the  main  result  of  this  part  of  the  project.  It  applies  to  any  two  -phase  con¬ 
servative  medium  with  a  well-defined  Lagrangian  density.  We  have  also  generalized  it  to  include 
dissipation  by  using  a  generalization  of  Hamilton’s  principle.  The  derivation  and  some  applica¬ 
tions  to  different  models  for  L  are  in  a  completed  paper  [1]  (see  Chapter  4  below).  Concerning 
the  specific  form  of  Eq.  (2.1)  we  note  that  the  left  hand  side  is  uncinticipated  emd  has  not  been 
previously  reported.  It  obviously  reduces  to  the  standard  mass  x acceleration  for  a  single  phase 
medium,  either  solid  or  fluid.  The  explicit  formula  for  the  interaction  force  is  entirely  new  and 
remains  to  be  properly  understood.  For  full  details  see  Chapter  4. 


2.3.3  Homogenization:  Precise  equations  for  layered  media 

As  a  first  step  in  understanding  the  MS  for  permeable  poroelastic  media  we  derived  the  equations 
governing  fimte  amplitude  motion  in  a  layered  fluid-solid  system.  The  layered  medium  is  the  only 
configuration  for  which  exact  averaging  yields  explicit  or  semi-explicit  analytical  expressions  for 
the  averaged  continuum.  The  stratified  medium  is  also  the  simplest  realization  of  a  permeable 
porous  medium,  and  has  been  previously  studied  both  theoretically  and  experimentally  (see  [2] 
for  references).  It  is  rich  system  for  linear  wave  motion  because  it  displays  both  the  Biot  fast 
and  slow  waves,  and  they  are  both  anisotropic. 

We  successfully  applied  the  method  of  homogenization  to  derive  the  MS  and  the  nonlinear 
constitutive  equations  starting  from  the  fundamental  equations  of  elasticity  and  fluid  mechanics 
for  the  constituents.  We  find  that  the  averaged  continuum  is  a  transversely  isotropic  Biot-type 
medium  with  a  fully  nonlinear  elastic  strain  energy.  The  full  details  eire  in  a  paper  [2]  which, 
for  the  first  time,  presents  explicit  nonlinear  Biot-like  equations  for  a  specific  microgeometry. 
It  also  gives  an  exact  form  of  the  permeability  convolution  operator  which  reflects  the  viscous 
effects  of  the  pore  fluid.  Previous  studies  have  derived  this  in  the  frequency  domain,  which  is 
not  sufficient  for  application  to  nonlinear  problems.  Our  time  dependent  permeability  operator 
displays  a  new  type  of  singularity  (~  1/v^)  which  must  be  included  in  shock  wave  models. 
We  also  examined  in  depth  the  linearized  theory.  Despite  several  previous  studies,  we  are  the 
first  to  present  the  complete  set  of  lineiirized  Biot  equations  for  the  stratified  fluid/solid  system. 
The  linearized  equations  possess  fast  emd  slow  wave  solutions,  which  have  been  compared  with 
ultrasonic  experiments  performed  in  the  1980’s.  This  work  has  been  published  [2],  and  a  copy  is 
attached  to  the  Report. 


2.3.4  Acoustoelasticity  for  fluid-saturated  porous  media 

One  of  the  simplest  nonlinear  phenomena  is  the  elastoacoustic  effect  whereby  the  speeds  of  small 
amplitude  waves  are  changed  by  applying  stress  or  strain.  The  theory  of  acoustoelasticity  has 
been  thoroughly  discussed  for  purely  elastic  materials  said  is  now  being  used  in  geoacoustics  to 
quantify  the  nonlineeirity  of  rocks  and  granular  media  (see  e.g.  [5]).  As  part  of  this  project  we 
derived  the  general  theory  for  small  dynamic  motion  superimposed  upon  large  static  deformation 


CHAPTER  2.  SUMMARY  OF  RESEARCH  ACCOMPLISHMENTS 


10 


in  isotropic  fluid-filled  poroelastic  solids.  Formulae  were  obtained  for  the  change  in  acoustic  wave 
speeds  for  arbitrary  loading,  both  on  the  frame  and  the  pore  fluid. 

Based  on  the  equations  we  obtained  for  the  change  of  the  fast,  slow  and  shear  wave  speeds  we 
proposed  a  set  of  specific  experiments  to  find  the  complete  set  of  third  order  elastic  moduli  for 
an  isotropic  poroelastic  medium.  Because  of  the  larger  munber  of  third  order  moduli  involved,  7 
as  compared  with  3  for  a  simple  elastic  medium,  experiments  combining  open-pore,  closed-pore, 
jacketed,  cind  unjacketed  configurations  are  required.  The  details  for  each  type  of  loading  are 
presented  and  a  set  of  possible  experiments  is  discussed.  FuU  details  can  be  foimd  in  [3]  (see 
Chapter  5  below). 


2.3.5  Waves  and  flutter  in  fluid  conveying  elastic  media 

This  part  of  the  work  was  undertaken  to  understand  the  stability  of  small  dynamic  motions  in 
porous  media  when  there  is  a  flow  present  in  the  channels.  As  a  simple  but  illustrative  case 
study,  we  considered  a  fluid  flowing  through  a  periodically  layered  system  comprising  alternating 
elastic  solid  and  liquid  constituents.  An  exact  analysis  was  performed  to  consider  the  existence 
and  stability  of  small  acoustic  waves  and  disturbances.  The  analysis  indicates  that  the  presence 
of  the  flow  introduces  the  possibility  of  flow  induced  flutter.  Unstable  waves  are  generedly  possible 
for  M  of  order  unity,  where  M  is  the  she2u:  wave  Mach  number,  but  appear  for  much  lower  values 
for  asymmetric  flexural  type  motion.  In  that  case  it  is  foimd  that  a  critical  wavenumber  exists, 
indicating  that  the  system  is  inherently  unstable  to  long  wavelength  disturbances.  For  full  details 
see  Chapter  6  below. 


2.4  Current  and  Future  Work 

Now  that  we  have  established  the  general  form  of  the  MS  in  Eulerian  variables  [1]  (Chapter  4 
below)  and  have  some  sense  of  the  form  of  the  nonlinear  elastic  strain  energy  and  permeability 
operators  [2]  we  believe  the  work  described  in  Section  2.3.1  is  the  most  importctnt.  Further 
numerical  experiments  will  result  in  better  insight  of  the  properties  of  the  MS  and  shock  waves 
in  media  having  both  Darcy-like  and  Navier-like  dissipation.  The  numeric-experimental  results 
will  show  us  the  essentially  novel  features  of  shock  wave  structure  in  fluid  permeable  elastic  media. 
Our  ultimate  objective  is  to  confirm  these  features  by  analytical  methods  as  much  as  possible. 
The  breakthrough  with  the  Hugoniot  conditions  will  open  wide  opportunities  in  understanding 
physiczil  processes  accompanying  nonlinear  wave  propagation  in  fluid  permeable  elastic  media  - 
we  believe  that  we  are  close  to  this  goal. 


Chapter  3 

in  two-phase  fluids 


finite  amplitude  waves  in  a  two  phase  liquid  medium 
with  viscous  damping.  The  equations  of  motion  are  considered  in  Eulerian  variables,  and  are 
based  upon  a  recently  derived  Hamiltonian  formulation  for  two  phase  media.  Two  distinct  types 
of  shocks  can  propagate,  analogous  to  fast  and  slow  waves.  The  speed  of  a  ID  shock  front  of 
either  type  is  determined  by  two  simtdtaneous  equations.  A  Hugoniot  condition  for  the  change 
#  in  the  internal  energy  is  derived.  A  numerical  exeimple  illustrates  the  shock  profile  for  the  two 

wave  modes. 


Shock  structure 

Traveling  wave  solutions  are  obtained  for 


3.1  Introduction 

One  might  suppose  that  the  appropriate  Euler-Lagrange  equations  must  be  well  known  for  a 
two  phase  medium.  But  the  answer  is  apparently  in  the  negative.  In  order  to  tmder stand  the 
problem,  let  us  first  consider  a  ID  single  phase  medium,  with  *  the  current  (Eulerian)  coordinate, 
u(x,t)  the  particle  displacement,  L{x,u,Ux,Ut)  the  Lagrangian  density  per  unit  volume.  The 
Euler-Lagrange  equation  is  then 

d  dL  dL 

dt  dut  ^  dx  dux  du  (^•^) 

Note  that  the  total  particle  velocity  is  w  =  ^  =  or  v  =  (1-Ua,)“^ut.  Also,  the  current 

mass  density  is  p  =  /!>o(l  -  «*)>  where  po  is  the  initial,  undeformed  density.  Let  E{ux)  be  the 
stored  energy  per  unit  mass,  then  the  Lagrangian  density  per  unit  volume  follows  as 

1  2 

I  =  (3.2) 

Substitution  of  (3.2)  into  (3.1),  yields  after  some  simplification,  the  well-known  equilibrium 
equation 


Dv 


P 


Dt 


+ 


(3.3) 
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where  the  pressure  p  is  defined  as 

2dE{p) 

P  -  P  —3 — . 
dp 

and  E{p)  =  E{ux).  The  point  here  is  that  simple  equations  are  not  simple  to  derive  using 
Euler-Lagr2mge  equations  in  current  coordinates. 

Another  way  of  phrasing  the  issue  is  that  the  fined  equation,  (3.4),  involves  quemtities  such 
as  the  toted  velocity  v  and  the  current  density  p  as  the  natural  pareuneters,  whereas  the  Eider- 
Lagreinge  equations  eire  defined  by  the  more  fundemiental  peirameters:  u  and  its  partial  derivatives 
with  respect  to  space  and  time.  In  fact,  the  E-L  equations  eire  derived  by  taking  independent 
Veiriations  of  these  quantities.  More  generally,  the  natural  p2irameters  of  study  in  field  theories 
of  continuum  mechanics  are  quantities  such  as  the  velocity  v  =  uj  +  v.Vu,  the  density  p  = 
p{det{l  -  Vu)),  and  the  deformation  gradient  Vu  in  elasticity  theories^  The  crucial  difficulty  is 
that  V,  p  2ind  Vu  are  not  independent,  and  it  is  not  immediately  clear  how  they  may  serve  as 
the  basis  for  a  variational  formulation.  That  is,  the  statement  of  the  E-L  equations  in  terms  of 
these  peirameters  is  non-trivial,  but  requires  tedious  algebra,  even  though  the  Lagrangian  L  is 
naturedly  defined  by  these  variables:  L  =  \pv^  —  pE{p)- 

The  situation  is  edtogether  of  a  different  order  of  difficulty  for  two-  and  multi-phase  media. 
For  instance,  consider  a  theory  for  a  fluid- saturated  porous  medium  in  which  we  posit  or  assume 
two  sets  of  physical  variables  at  the  point  x,  see  Table  1. 


Table  1 


Medium 

Mass  Density 

Velocity 

Strain 

Solid 

Ps 

Va 

Vu, 

Fluid 

pf 

V/ 

Vu/ 

Table  1.  The  physical  parameters  required  for  a  description  of  a  solid/fiuid  continuum.  The 
subscripts  s  and  f  refer  to  the  solid  emd  fluid  phases,  respectively. 

The  most  genercd  form  of  the  Lagrangi2in  density  per  unit  total  mass  is  therefore  of  the  form 

L  =  L  (p„  p/,  v„  V/,  Vu„  Vu/)  .  (3.5) 

^Note  that  the  deformation  gradient  as  defined  here  is  not  the  same  as  that  normally  encountered  in  nonlinear 
solid  mechanics,  but  they  may  be  easily  related  to  one  another. 
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How  is  one  to  obtain  the  equations  of  motion  from  this?  The  answer  is  not  at  all  obvious,  because 
the  petreuneters  vy,  Vu,,  and  Vu/  are  not  independent  quantities.  They  are  related 

by  their  defimtions  cind  by  the  physical  constraints  of  mass  conservation  for  each  phase,  which 
imply  the  pointwise  identities 

^  +  V.(p„v„)  =  0  (3.6) 

Recently,  we  have  shown  how  to  derive  the  governing  equations  by  using  Hamilton’s  principle 
as  the  stcirting  point  [1].  The  derivation  is  detailed  and  imedifying,  but  relies  upon  the  functional 
relations  between  Pa»Va>  Vua  to  derive  a  lineeir  form  in  the  virtual  velocities  f)j  =  + 

f^.Vua,  where  t  is  the  parameter  of  variation  in  the  Principle  of  Least  Action.  This  gives  two 
equations:  one  for  the  solid  and  one  for  the  fluid,  each  of  the  form 

~  P  ^Pb^iPb  +  ViVjUfc5 j 

+V,(/)p6lrpJ  -  ^jipLvju^^A'kib),  b  =  SOT  f,  (3.7) 

where 

^=f  +  v,.v  4,.(x.()  =  «;-VX,  (3.8) 

and  p  =  p,  pf  is  the  total  mass  density  of  the  system. 

Heunilton’s  Principle  of  Least  Action  may  he  extended  to  account  for  irreversible  effects, 
specificeJly,  viscous  body  forces  and  viscous  stresses.  We  also  considered  such  irreversible  effects 
cind  derived  the  governing  equations,  which  turn  out  to  be  similar  to  (3.7)  except  that  the  right 
members  are  supplemented  by  the  terms 

^ jPLi,  +  Pb.vi»  ^  b  =  soTf,  (3.9) 

where  viscous  forces,  and  viscous  stresses. 

In  summary,  these  governing  equations  are  fully  general,  and  are  suitable  for  nonlinear  dy¬ 
namic  modeling  of  two  phase  media.  The  equations  depend  upon  the  explicit  form  of  L  used, 
but  remain  too  general  and  too  complex  for  studying  shock  structure  in  fluid/solid  continua. 
Instead,  we  consider  here  the  simpler  two  phase  medium:  a  one  dimensional  model  of  a  two 
liquid  medium. 

The  two  liquid  medirun  still  displays  many  of  the  most  relevant  characteristics  of  fluid- 
saturated  porous  solids,  such  as  the  existence  of  two  distinct  compressional  wave  types  for 
infinitesimal  motions.  The  structure  of  shocks  waves  in  uniform  fluids  is  well  understood  by 
now,  e.g.  [49],  but  there  is  no  corresponding  theory  available  for  systems  of  two  or  more  phases 
(models  exist  which  include  the  influence  of  subsidiary  phases,  such  as  secondary  gases  in  gas 
mixtures,  but  they  are,  as  far  as  we  understand  them,  essentially  single  phase  theories).  Among 
our  objectives  here  are  the  following:  (i)  to  determine  if  two  types  of  shocks  are  possible  (fast  and 
slow),  (ii)  to  tmderstztnd  the  dependence  of  shock  speed  on  amplitude,  and  (iii)  to  understand 
the  role  of  viscous  effects.  The  permeability  is  another  important  parameter  in  the  dynamics  of 
porous  media,  and  cm  be  included  here  as  a  dissipative  (viscous)  term.  However,  only  ‘direct’ 
viscosity,  will  be  considered  here.  Also,  we  ignore  entropy  effects  in  this  paper. 
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3.2  Governing  ID  equations  of  a  liquid  mixture 


The  physic2tl  quantities  of  interest  cire  the  densities  pi,  p2y  and  the  particle  velocities  Vi^  V2,  all 
of  which  are  functions  of  time  and  the  single  spatial  coordinate  x.  The  total  mass  density  at  a 
point  is  p  =  Pi  +  P2-  The  equations  of  mass  conservation  for  each  species  axe 


DiPi  dvi 


D2P2  ,  dv2 
Dt  ~  ’ 


where  the  total  material  derivatives  for  each  phase  are 


Di  d  d 
Dt  dt 


D2  d  ,  d 
Dt  dt  dx 


(3.10) 


(3.11) 


The  governing  equations  of  motion  may  be  expressed  succinctly  in  terms  of  Lagrangian  density 
per  unit  total  mass  L{vi,V2,Pi,P2),  or  explicitly 


L  =  T{vi,V2,PuP2)  -  E{pi,p2), 


(3.12) 


where  E  is  the  elastic  energy  function,  and  T  the  kinetic  energy.  Thus,  from  eqs.  (3.7)-(3.9)  and 

[1], 


(  P  rp  \ 

^^Dt 

= 

D2  , 

(  P  rr  \ 

^^Dt  * 

U2  “ 

where  the  partial  pressures  are 

dpi  _ 

-■^  +  ^12  +  ^1-^, 
4.  F  4.  ,, 


Pi  =  -/>Plbp,,  P2  =  -PP2bp2, 

and  Fi2,  F21  may  be  interpreted  as  equal  but  opposite  interaction  force  densities, 
Fi2  = 

The  total  mechanical  energy  density  (elastic  plus  kinetic)  is  [1] 

S  —  “1“  ^2^V2  > 

and  the  toted  energy  flux  at  a  point  is  [1] 

T  -PlVi-^  P2V2  +  (Ui  -  V2)p2VlLy^  +  {V2  -  Vi)piV2L^,^  . 


dv2 

dx 


dp2 

dx 


F21  =  PiLv2  +  PlI'P2  “  P2b„,  -  P2bp, 


dvi 

dx 


dpi 

dx 


(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 


The  equation  for  the  balance  of  energy  at  a  point  follows  from  the  equations  of  motion  as  [1] 


DiS  D2S  ,  dF  fdviY  (dv2\^ 

'’■‘m ■ 


(3.19) 
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3.3  Propagating  waves 

3.3.1  Shock  speed  and  modes 

Consider  a  propagating  wave  of  the  form  f(x,t)  =  f(^)  where  ^  =  a:  -  ct  for  any  veiriable  f. 
Thus,  Daf/Dt  =  (va  -  c)f',  and  so  the  mass  conservation  equations  become 

{c-vi)pi  =  ai,  (c  -  V2)/>2  =  02,  (3.20) 

where  oi  and  02  eire  constants.  Using  (3.20)  to  simplify  F12  and  F21,  the  equations  of  motion 
(3.13)  become 


-Pi'  +  Pi^-P2^  +  Aii<  =  0, 

PPl  PP2 

(3.21) 

//  ^  Q 

PPl  PP2 

(3.22) 

where 

Pa  —  Pa  ^a  Pva  • 

Pa 

(3.23) 

Note  that 

Pa  —  PPaPpa  ) 

(3.24) 

where 

L{pi,p2)  =  l(c-^,C-^,Pl,P2]  , 

\  Pi  P2  J 

(3.25) 

is  a  function  only  of  pi  and  p2  because  vi  and  V2  have  been  eliminated  using  (3.20). 
Adding  the  two  equations  of  motion  and  integrating,  implies  that 


+  P2V2  -  Pi  -  P2  =  b,  (3.26) 

where  b  is  another  constant.  The  constants  can  be  determined  from  the  conditions  at  ^  =  ±oo. 
The  wave  is  assumed  to  propagate  into  an  undisturbed  medium,  implying  0,  as  f  oo, 

and  hence 

~  Ploo^j  0>2  “  P2oo^y  b  =  “Ploo  ^  P2ooy  (3.27) 

where  piooj  P200  aJid  pioo?  P200?  are  the  equilibrium  densities  and  partial  pressures  ahead  of  the 
wave.  NormaJly,  and  in  the  following,  both  pioo  and  p2oo  are  zero.  The  densities  at  2iny  point 
can  now  be  expressed  in  terms  of  the  velocities  using  equations  (3.20)  and  (3.27), 


Ploo^ 

Pi  = 

P2oo^ 

P2  =  , 

(3.28) 

C  —  Vi 

C  —  Vi 

~  (Pl  *“  Ploo)  ) 

V2  =  {P2  -  P200)  —  ■ 

(3.29) 

Pi 

P2 

Pa  —  Va  *“  Poo  C  Pva  > 


Also,  note  that 


(3.30) 
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where  poo  =  Pioo  +  P2oo-  Integrating  the  fundamental  equations  (3.21)  yields 

Plivi.gh,V2.sh)  =  Ploo  +  IW  ,  P2{V\.»hjV2.sh)  =  P2oo  -  Iw  ,  (3.31) 


where 

OO  ^  / 

Iw  =  -  1  Fndi  =  /  di 

J  J  \  PP2  PPl  ) 

— OO  — OO 

OO 

=  J  {p2lpiPi-  PiIp2P'^  d^.  (3.32) 

—  OO 

The  pair  of  equations  (3.31)  provides  a  nonlinear  eigenvalue  system  for  the  wave  speed  c  and  the 
pmticle  velocities  v\,sh  and  V2.,h-  The  system  (3.31)  as  it  stands  is  not  closed  but  requires  prior 
knowledge  of  the  integreil  of  the  interaction  force  Iw- 

3.3.2  A  Hugoniot  relation 

When  applied  to  traveling  waves  the  energy  balaince  (3.19)  becomes,  using  the  identities  (3.29), 

-p^c  -Ml  -  P2  .  (3.33) 

Integrating  through  the  shock  yields 

Poo^  ~~  Poo^^oo  —  "h  ^vig  ?  (3.34) 

where  Dyis  is  the  total  viscous  dissipation  for  the  shock, 

OO 

D,i,  =  J  [mi  {v[f  +  M2  (^)']  d^  -  (3.35) 

—  OO 

The  elastic  energy  density  2Lhead  of  the  shock  is  assumed  to  be  zero  with  no  loss  in  generality, 
5oo  =  0.  Also,  the  flux  far  behind  the  shock,  !Fthi  can  be  simplified  using  the  characteristic 
equations  (3.31).  Thus,  (3.34)  becomes 

MooC  Dyxs  —  P\<x>^'^\,ghTvi  "i"  P2oo^'^2,»hTxj2 


+(vi  -  V2)  {p2.shVi.3hTvi  -  pi.»hV2.»hTv2  +  Iw)  -  (3.36) 

The  left  member  reflects  the  total  change  in  internal  energy  density  caused  by  the  passage  of  the 
shock.  This  Huguniot-type  relation  also  connects  the  integral  Iw  with  the  viscous  dissipation 
integral 


CHAPTER  3.  SHOCK  STRUCTURE  IN  TWO-PHASE  FLUIDS 


17 


3.4  Shock  structure 

The  structure  of  the  shock  follows  from  the  ODEs  (3.21).  Alternatively,  taking  a  lineair  super¬ 
position  of  the  two  equations  of  motion  (3.21)  gives 

PiM  -  v'')  +  (Ml  +  M2)  (Pi—  -  P2^)  +  P1P2  -  f^iPi  =  0 .  (3.37) 

\  PPl  PP2  ) 

We  will  use  the  alternative  pedr  of  ODEs  (3.26)  and  (3.37)  which  define  a  third  order  system  as 
compared  with  the  original  fourth  order  equations  (3.21).  Define 

U  =  Mi^i  +  P2V2  y  (3.38) 

V  =  V1-V2  (3.39) 

w  '  '  P^  ,  P2 

Pi  P2 


then  we  may  rewrite  equation  (3.37)  as 


=  (fh'  +  /•?')  -  Pi^)  ■ 

V  PP2  PPl  / 


(3.40) 


(3.41) 


The  derivatives  pi  aind  P2  can  be  expressed  in  terms  of  derivatives  of  U  ^lnd  V  by  using  (3.28) 
zmd  the  relations  vi  =  (17  -f  M2V^)/(Mi  +  ^2)*  ^2  =  {U  —  MiV')/(mi  +  M2)-  The  equations  therefore 
reduce  to  a  system  of  three  first  order  ODEs: 


U'{i)  =  El  -|-  P2  -  Ploo  -  J>2oo  > 


(3.42) 


no  =  — +w, 

Pi  P2 


(3.43) 


W'(0  =  -^^\—{U'-firV')-^{U'  +  H2V')\. 

/il/X2/>C  [p2oc  Ploo  J 

Behind  the  wave  the  velocities  attain  non-zero  asymptotic  values,  Thus,  Va 
-00.  The  densities  far  behind  the  wave  follow  from  (3.28)  as 


Puh  = 


c  -  Vi,sh 


P2$h  = 


C  -  V2.sh 


Equations  (3.42)  and  (3.43)  therefore  imply  the  identities 

Pl{vi.3h,V2,sh)  +  P2{'^l,»h,V2.sh)  =  Ploo  +  P2oo  ? 


(3.44) 


Va.sh  as  ^ 


(3.45) 


(3.46) 


Pli,'^l.shy  ^2.g/t)  _  ^2(^1.*/^  '^2.j/t)  _  Ploo  _  P2oo  .  U™  pp-  _  pp 

Ml  P2  Pi  P2  i-*-oo 


(3.47) 
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This  pair  of  equations  Ccin  be  rewritten  as  (3.31)  because  of  the  identity 

OO 

Iw  =  (/ii"'  (3.48) 

—  OO 


In  the  linearized  limit,  the  terms  pi 
eigenvalue  system  reduces  to 

'^2,sh) 


OO,  p2oo  W  are  second  order,  and  therefore  the 
=  0,  P2(t'i.*fc,V2.»/.)  =  0.  (3.49) 


Explicitly  linearizing  yields  a  pair  of  simultaneous  equations 


Pl^PlPl 

PlP2Ep^P2  C^Ty^y^ 

pVl.sh 

■  0  ■ 

PlPlEp^p^  —  C^Ty^y^ 

P2EP2P2  " 

pn.th  _ 

_  0  _ 

(3.50) 


3.4.1  Single  phase  medium 

For  a  single  phase  medium,  equation  (3.42)  becomes  “Burger’s”  equations,  that  is, 


pv'  =  -Poo  -  P^Lp , 


or  in  terms  of  p  alone. 


p'  .  .2  Poo) 

PPoo^~2  —  P\P)  ~  Poo  ~  Poo^ 

P  P 


Equation  (3.52)  can  be  integrated  as 


^  =  P  I  \—{Poo-p)  +  :^iVoo-v) 
J  I  C  Tfoo 

»?(0) 


-1 


dr), 


(3.51) 

(3.52) 


(3.53) 


where  7}  =  1/ p. 

The  characteristic  equation  is  the  single  condition  (3.46),  which  becomes  (see  p.  159  of  [50]) 


2  P»h  ( Psh  Poo\  2  _  Poo){psh  Poo) 

C  -  I  -  I  J  =  - . 

Poo  \Psh  “  PooJ  PshPoo 

Therefore,  eliminating  c  from  (3.53)  yields 

_  f  _ \/PshPooiPsh  Poo)iPeh  ~  Poo)  ^P _ 

P  J  pHp{p)- Poo)iP»h- Poo)  -  PPsh{P,h-Poo){p- Poo) 
P(0) 


(3.54) 


(3.55) 


The  constant  of  integration  can  be  eliminated  by  defining  ^  =  0  to  be  the  point  at  which  p  has 
a  certain  value  between  poo  and  p,h-  For  instance,  we  may  set  p(0)  =  j(pj,h  +  ^oo),  then  (3.55) 
yields  the  shock  profile. 
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As  a  first  approximation  to  the  exact  shock  structure,  the  denominator  in  the  integral  can 
be  replaced  by  a  polynomial  with  zeros  as  p  =  p^h  and  p  =  poo-  The  simplest  approximation  is 
a  quadratic, 

p^ipip)  -  Poo){Psh  -  Poo)  -  pp,h{p»h  -  Poo){p  -  Poo)  «  p,h{p,h  -  Poo){p  -  p.h){p  -  Poo) ,  (3.56) 
which  can  be  integrated  exactly,  to  give 

~  ~  tanh  —  ^^^2^  (3.57) 

3*4.2  Linearized  equations 

In  the  linear  approximation  the  terms  involving  F12  and  F21  of  (3.16)  are  of  second  order  in  the 
small  quantities:  the  deviations  in  the  densities  and  the  particle  velocities.  Let  vi  and  V2  be  the 
independent  variables,  then  the  linearized  equations  are,  from  (3.21), 

//iVi  +  pVi  +  ^cTvit;2  ~^PlP2^  P'^2 

p2'^2  “t“  (^'^ViV2  “  “®P1P2^  P^l  “I"  ^-®P2P2^  P^2 

3*5  Numerical  experiments 

The  numerical  problem  can  be  cast  as  three  coupled  first  order  ODE’s  for  Vi  and  vj,  defined 
by  eq.  (3.26)  with  6=0  and  eq.  (3.37).  This  presents  a  nonlinear  eigenvalue  problem  for  c, 
and  the  values  of  Vi  and  V2  (or  pi  Euid  P2)  behind  the  shock  cue  Vi  V2  V2.sh’,  as 

^  —00.  The  nonlinear  eigenvalue  problem  can  be  expressed  in  almost  cJgebraic  form  via  eqs. 

(3.31)  with  pioo  =  0  and  p2oo  =  0.  The  difficulty  is  due  to  the  term  Iw,  which  is  the  total 
integrated  interaction  force  through  the  shock.  For  single  phase:  7^^  =  0  and  we  immediately 
have  an  algebredc  equation  for  c  =  c(t;,h)*  For  two  phases:  c  =  c(vi.,/i,  t;2.5h»-ftv)>  which  requires 
explicit  numerical  solution.  For  infinitesimal  motion,  we  have  Iw  =  0,  and  thus  obtain  two  linear 
equations  for  V2.Bh  (Biot-type  equations)  with  fast  and  slow  roots  for  c. 

Numerical  experiments  were  performed  to  solve  for  the  shock  speed  and  structure.  The 
procedure  used  (so  far)  is  to  solve  the  ODES  iteratively.  First,  fix  vi(0),  guess  ^2(6)  cind  c.  Then 
integrate  from  ^  =  0  both  ways,  and  subsequently  adjust  1^2(0)  and  c  to  iteratively  reach  steady 
states  at  large  z.  We  find  that  there  can  exist  two  types  of  shock  fronts:  fast  and  slow,  with  the 
following  chciracteristics.  the  fast  shock  has  >  0  (in  phase),  while  the  slow  shock  has 

<  0  (out  of  phase).  A  typical  pair  of  solutions  is  shown  in  Figure  3.1. 


0, 

(3.58) 

0. 

(3.59) 

Chapter  4 

Hamiltonian  methods  for 
fluid-saturated  porous  media 


A  new  system  of  equations  governing  nonlinear  dynamics  of  fluid  permeable  poro-elastic  media  is 
derived  on  the  basis  of  Hcimilton’s  principle  for  reversible  effects  and  the  Onsager-Sedov  approach 
for  irreversible  effects.  The  equations  Me  in  Eulerian  variables,  suitable  for  dealing  with  fluid  Euid 
sohd  phenomena  simultzineously.  The  classical  (Murnaghan-Hke)  equations  of  nonlineM  elasticity, 
as  well  as  the  governing  equations  of  the  ideal  cind  Navier-Stokes  fluids  Me  shown  to  be  special 
cases  of  the  general  governing  system.  It  is  weU-known  that  the  Navier-Stokes  equations  of 
compressible  fluid  provide  a  correct  self-consistent  basis  for  studying  a  wide  VMiety  of  nonhneM 
effects  in  fluids.  The  authors  believe  that  the  governing  system  proposed  here  provides  the  same 
opportunities  for  VMious  nonlineM  effects  in  poro-elastic  fluid-penetrable  media. 


4.1  Introduction 

It  is  common  practice  to  formulate  nonlineM  equations  of  solids  in  reference  or  Lagr2ingian 
coordinates.  Fluid  dyneimics,  on  the  other  hand,  is  more  naturally  considered  in  current  or 
Etderiein  coordinates.  The  reasons  Me  cleM:  in  solids  the  material  pMticles  adhere  to  one 
another,  so  that  the  reference  description  always  provides  a  continuous  mapping  to  the  current 
configuration  of  pMticles.  Not  so  for  fluids,  where  iMge  sccde  sheMing  crm  sepMate  neighboring 
pMticles  in  the  reference  system.  The  same  is  true  of  fluid-solid  composite  media,  because  the 
fluid  can  shde  past  sohd  pMticles.  The  use  of  Eulerian  coordinates  for  such  media  therefore 
appeMS  to  be  the  most  convenient  procedure,  and  maybe  the  only  feasible  one. 

The  purpose  of  this  Chapter  is  to  estabhsh  the  governing  equations  for  fluid-sohd  media  in 
Eulerian  coordinates  by  means  of  Hamilton’s  principle  for  conservative  media  zmd  the  Onsager- 
Sedov  technique  for  irreversible  effects.  The  idea  is  to  stMt  with  a  well  defined  Lagremgian 
density,  from  which  all  equations  should  follow.  Although  this  approach  is  potentially  the  most 
rehable  and  flexible  its  usage  for  the  Euleriem  description  of  continuous  media  demands  some  CMe. 
However,  we  beheve  this  is  the  most  general  and  consistent  approach  to  formulating  nonhneM 
equations  of  fluid-sohd  continua.  In  pMticulM,  it  allows  one  to  study  the  internal  structure  of 
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shock  waves  and  flutter-like  phenomena,  among  others. 

We  note  that  there  are  many  excellent  articles  on  the  fundamental  equations  of  fluid-permeable 
solids.  Nonlinear  treatments  are  usually  framed  in  terms  of  a  mixture  theory  [41,  42]  for  which 
the  starting  point  is  a  set  of  balance  laws  for  mass,  momentum,  etc.  Here  we  take  one  step  back 
and  invoke  a  stationcirity  principle,  that  of  least  action,  to  generate  the  eqtdlibrium  equations. 
Conservation  of  mass  is  considered  fundamental  for  each  species,  but  no  other  balance  laws  are 
invoked.  Mixture  theories  can  also  be  formulated  from  an  energy  stand  point  [43],  and  nonlinear 
poroelasticity  has  been  developed  from  a  thermodynamic  point  of  view  [44].  However,  we  are 
unaware  of  a  treatment  based  on  the  Hamiltonian  principle  of  least  action  as  developed  here. 

In  many  cases  discrete  and  distributed  mechanical  systems  show  a  behavior  which  can  be 
treated  as  conservative.  It  seems  reasonable  to  study  such  behavior  using  a  proper  generalization 
of  the  Hamilton  variational  principle.  Following  some  preliminary  deflnitions  in  Section  2  we 
begin  in  Section  3  by  considering  the  procedure  of  [45]- [47]  for  performing  variations  in  Eulerian 
vetriables.  The  remainder  of  Section  3  describes  the  derivation  of  the  governing  equations  for 
a  conservative  system  from  Hamilton’s  principle.  We  build  upon  these  results  in  Section  4  and 
include  dissipative  or  irreversible  effects,  associated  with  fluid  viscosity  for  instance.  The  starting 
point  is  a  generalization  of  Hamilton’s  principle  using  the  Onsager-Sedov  method.  Some  general 
results  concerning  energy  are  discussed.  Finally  in  Section  5  we  provide  a  variational  derivation 
of  the  dynamic  equations  of  classiccJ  hydrodynamics  and  of  Murnaghan’s  equations  of  nonlinear 
elasticity  [48].  We  also  present  three  examples  of  the  application  of  the  governing  equations  to 
two  component  media  included  fluid-flUed  poroelastic  solids. 

4.2  Notation  and  the  Eulerian  description  of  two-component 
continua 

The  use  of  the  Euleristn  description  has  several  advantages  for  analyzing  two-component  media, 
and  fluid  permeable  solids  in  particular.  In  the  Eulerian  description  time  t  and  the  coordinates 
2*  (Latin  indices  like  i^j^kj  take  the  values  1,2,3)  are  used  as  independent  Vciriables.  For  the 
sake  of  simplicity  the  reader  Ccin  think  that  the  system  is  referred  to  the  Cartesian  coordinates 
although  we  prefer  to  use  a  covariant  form  of  the  equations  (in  particular,  this  form  makes 
general  equations  much  more  eloquent  and  expressive,  and  it  allows  one  to  switch  easily  between 
geometrically  different  coordinate  systems).  We  use  the  notation  for  the  co-  and 

contra- variant  metrics  and  covariant  differentiation  (in  the  Cartesian  system  Zij  =  z^^  =  Sj  and 
CO-  and  contra-variant  components  appezir  to  be  equal  while  coveiriant  differentiation  reduces 
to  the  partial  Vj  =  d/dz*).  In  order  to  distinguish  between  different  characteristics  of  the  two 
components  we  use  extra-subscripts  a,  6  taking  Vcdues  1  or  2  or  the  subscripts  s  or  f  when 
emphasizing  the  solid  or  fluid  nature  of  the  component.  In  particular,  we  use  the  notation 
u^(2,t),  pa{z^t)y  and  t;*(2,t)  for  the  displacements,  densities  and  velocities  of  finite  deformations 
of  the  components.  Introducing  the  “material”  time  derivative  Da/Dt  =  d/dt  +  v^^Vi  associated 
with  the  “a-th”  constituent  we  get  by  definition 

Daui 


d 
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The  following  pair  of  identities  zire  inunediate  consequences  of  Eq.  (4.1) 


dt 


dui{z,t) 

dt 


B\ 


(4.2) 


where  A\j^[z,t)  =  -  Vju’  and  B]j^{z,t)  is  the  inverse  of  A'^^,  i.e.,  B'^  A^/^  =  S^. 

The  fields  pa{z,t)  eind  v^{z,t)  obey  the  well-known  formula  of  mass  conservation  for  each 
phase, 

^  +  ^i{Pavi)  =  0.  (4.3) 

The  difFerenticil  form  of  the  mass  conservation  result  (4.3)  is  equivalent  to  the  following  algebraic 
one 


Pa  =  pldet\A'jJ, 


(4.4) 


where  is  the  mass  density  in  the  reference  configuration.  Differentiating  Eq.  (4.4)  we  get  the 
following  relationship 


dpa 

dViuL 


(4.5) 


4.3  The  Hamiltonian  approach  for  two-component  media  in  the 
Eulerian  variables 

Assuming  the  absence  of  externcil  mass  and  surface  forces  we  choose  the  Action  S  corresponding 
to  a  “cold”  multi-component  substance  in  the  following  standard  form: 

S-  f  dt  f  dQ.p{zjt)L^  (4.6) 

JtQ  Jil 

where  L  =  T  -  E  and  T  and  E  are  the  kinetic  and  internal  (elastic)  energies  of  the  substcince 
per  unit  mass,  and  p{z,t)  is  the  mass  density  of  the  substance  which  is  equal  to  the  sum  of  the 
partial  densities  pa{z^t)  of  the  constituents 

a 

In  what  follows  we  assume  that  the  kinetic  energy  density  is  an  algebraic  function  of  velocities 
of  the  constituents  t;*(z,f)  while  the  internal  energy  density  is  a  function  of  the  displacement 
gradients  VjU^(z, t)  and  the  partial  densities  pa{Zjt): 

T  =  T{vi),  E  =  £;(VX,  Pa),  L  =  L{pa,  < ,  V.V)  •  (4-8) 

The  present  2Lncilysis  applies  to  ^lny  number  of  constituents,  a  =  1, 2,  •  •  •,  but  for  practiced  purposes 
we  win  be  concerned  with  two- component  media. 

Dyneunic  equations  and  natural  boimdary  conditions  of  the  actual  motion  cein  be  obtained  as 
the  conditions  of  veinishing  of  the  first  variation  of  the  Action  S  for  aU  “admissible”  trajectories 
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In  addition  to  some  natnrzd  demands  of  sufficient  smoothness,  aJl  ad- 
missible  fields  have  to  obey  certain  “mechanical’  constraints.  First  of  afi,  as  always,  aU  admissible 
dynamical  fields  have  to  coincide  with  the  actual  fields  in  the  initial  (at  t  =  to)  and  final  (at 
t  =  ti)  configurations.  Other  crucial  constraints  follow  from  the  fact  that  the  fields 
are  mutually  dependent  since,  say,  the  two  former  fields  can  be  foimd  from  the  latter  by  means 
of  spatial  and  time  differentiation.  This  fact  implies  certain  linear  differential  relationships  for 
the  variations  of  these  fields,  and  the  independent  variations  should  be  carefully  extracted  in  the 
first  variation  of  the  Action  ^5.  In  order  to  make  this  point  more  transparent  let  us  consider  a 
one-parameter  family  of  the  admissible  fields 

Pa  =  Pa{z,  ty  t),  Vl  =  <(2,  t,  t),  u‘  =  (z,  t,  t),  (4.9) 


where  r  is  the  pareimeter  of  variation.  By  definition,  the  Euleri2m  variations  of  these  fields  are 
the  following  functions  of  z  and  t: 


Spaizyt)  =  ^iz,t,T) 


T=0 


—  Ra{Zy  t). 


6vi{z,t)  = 


=  Giiz,t), 


(4.10) 


It=0 


Su\{z,t)  =  ^(z,t,T) 


=  Qa{z,t). 


lr=0 


With  the  help  of  Eqs.  (4.2),  (4.4)  and  (4.5)  one  can  easily  establish  the  following  explicit  formulae 
for  the  variations  Ra{z,t)^Ga{z,t)  in  terms  of  the  variation  Qa{z^t)  emd  its  derivatives: 


G'^{z,t)  =  B':, 


dQ^ 

dt 


dvP 


RaiZyt)  —  PaB^.ia^jQa- 


(4.11) 


In  order  to  establish  a  system  of  governing  equations  one  ought  to  substitute  the  one- 
dimensioned  fields  (4.9)  into  (4.6),  differentiate  the  integral,  eind  then,  extract  the  independent 
vetriation  Q\  using  the  relationships  (4.11)  and  integration  by  parts  (in  space  eind  time).  Then, 
equating  to  zero  the  coefficients  of  in  the  integrand  we  arrive  at  the  desired  master  system. 
However,  this  approach,  although  absolutely  correct  conceptually,  requires  very  laborious  com¬ 
putation  when  the  Eulerian  description  is  used.  In  what  follows  we  use  a  somewhat  different 
technique  which  seems  to  be  much  more  convenient  when  dealing  with  the  Eulerian  description. 
In  this  technique  the  virtuzd  velocities  /‘(z,f,T)  of  the  constituents  eure  treated  as  independent 
variations  rather  than  the  functions  Su\{z,t).  The  virtual  velocities  /*(z,t,T)  are  defined  im¬ 
plicitly  by  the  relationship 

/a(^.  '^)  = 

where  Da/Dr  is  the  material  derivative  associated  with  the  parameter  r:  Da/Dr  =  d/dr+f^Vi. 
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The  following  statement  is  crucied  for  overcoming  future  computational  obstacles. 
Lemma.  The  operators  DalDt  and  Dal  Dr  commute  on  u^(z,t,T) 


Da  Da 

Dt'  Dr 


<  = 


Dg  Dg 

’m’Dr 


Dg  Dg  ^  j  _ 

“  Dt 


Dgfg 

Dt 


^  =  0. 

Dt 


Proof.  Differentiating  Eq.  (4.2)  we  get 


_  Dg  ( p»  \  _ 


D 

Dt 


D^ 

Dt 


dt 


n*  j_  O’ 

^  la  + 


dvlg  DgB\ 


ISL 


dt  Dt 


(4.13) 


(4.14) 


One  cam  easily  derive  the  following  relationships: 

DgBL 


Dt 

D 

Dt 

Dg^ 

Dt  dt 


D  4” 

oi  Dm 
-^.ng^.la~ 


■PgVm«a 


Dt 


Dt 

d  Dgul  dfi 


-  _4".  T7  fi 

—  ^.jg^rnlai 


dt  Dt 


dt 


dt 


Inserting  these  in  the  RHS  of  (4.14)  we  cirrive  at  the  desired  restdt 


DgV'a  _  df^a 


+  (V„./i)  BJa 


^«'g  ^/.g  I  A 


Dt  dt  dt 

We  remind  the  reader  of  the  following  formulae  for  a  derivative  of  the  integral  over  a  domain 
occupied  by  a  moving  body  consisting  of  particles  of  the  same  material: 


[  dIlpaL=  f  dQpa^^,  ^  [  d0.paL=  f  dUp, 

dT  Jn  Ju  Dt  dt  Jn  Jn 


DgL 
‘  Dt  ■ 


Combining  (4.6)  and  (4.15)  we  get  the  following  formula 
dT 


DgL 

Dt  ■ 


The  fundamental  equations  of  the  system  follow  by  equating  the  latter  to  zero,  which 
now  proceed  to  do.  First,  we  note  the  identity 

Dg  _  D^  ( fk  ■tk\Tj 

-^-D^-^fb-fgWk, 
which  combined  with  (4.13)  yields  the  relationships: 


(4.15) 

(4.16) 
we  will 

(4.17) 


Dgvi 


b  _ 


Dbfi 


Dt 

DgPb 

Dt 


Ub  -  fg)^kvl 


Dt 

-PbVkft  -  ifb  -  fa)^kPb, 


(4.18) 


^VpU,6  =  ^,r6Vp/?  -  (/t*' - /^)VfcVpU,6. 
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Hence,  the  integrand  of  Eq.  (4.16)  becomes 

^f>fb  ( tk  tk 


EDaL 

a  6 


L,.i 


^0  \  Dt 


ifb  -  a'^kvi 


+-^VpU,t  (AqrbVpfh  -  {ft  -  /^)VfcVpU,(,^] 

The  individual  terms  cem  be  simplified,  yielding  sequentially 

nT  £^IL  —  \^n£.T 

vi  Dt  Dt 


-  Lp,  {pb^kft  +  ifb  -  ft)'^kPb) 

(4.19) 


a  b 


a  b  b 

^  ^  [^P  {^P^^p'^qh^qrbfb^  ~~  fb^P  ]  > 

a  b  b 


(4.20) 


a  b 

where  Zkb  are  arbitrary,  and 


^  a 


is  the  “mean”  virtual  velocity.  Combining  (4.19)  and  (4.20)  we  get  the  relationship 


+ftVk{pPbLp,)  -  pLv,u,,^k^pU,b{ft  -  F^)  -  pLp,VkPb{ft  -  F’^) 


ek  vik\ 


-ft'^p{pLv,u,,Aq,b)  +  Vp(pivp.,»A,p6/?)  -  ^kipPbLpjt)]  .  (4.21) 

Now,  let  us  insert  Eq.  (4.21)  in  (4.16)  and  eliminate  all  spatial  divergence  terms  -  viz,  the  final 
two  in  Eq.  (4.21)  (here  and  in  the  following  we  tacitly  assume  that  aU  integrals  over  boundaries 
Vcmish:  the  reader  can  easily  figure  out  several  conservative  boundary  conditions  of  this  sort), 
yielding 
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-pL^Vkviift  -  F*')  +  ft^kipPbL,,)  -  pL,,VkPb{ft  -  F") 

o 


^fb^piP^'^P^qb-^qrb)  pUgb^ p'^qbifb  ^  )]  • 

Using  the  second  of  the  two  relationships  (4.15)  we  get 


t-h 


t=ti 


(4.22) 


(4.23) 


The  last  substitution  vanishes  because  the  initicd  and  final  configurations  of  aU  admissible  tra¬ 
jectories  are  the  same.  In  view  of  (4.23)  the  first  integral  in  Eq.  (4.22)  disappears.  Now,  by 
separating  independent  virtual  velocities  fl  in  (4.22)  and  setting  dS/dr  —  0,  we  obtciin  the 
following  governing  system  of  evolution: 

P^’^  ”  Pb^2  ^Pa^iPa^  ^  P^v^^^b  ""  P^Pb^iP^  pI^V ju^b^ j'^kb 

+  Vt(pP6ipJ  -  Vj(pi/V,u*6^ikt6)  +  PbY^L^jUko^i^j'^ka, 


which  we  can  rewrite  as 

(^'^’'0  ~  ^  (^vi^'^b  +  ^Pi^ipb  +  Lvjuti^i^jUkb) 

+Vi(/>/)bFpJ  -  jipLvjuttAkib)-  (4-24) 

This  is  the  first  significant  result  of  the  Chapter.  In  the  next  Section  we  will  generahze  this  to 
include  dissipation. 


4.4  The  Hamilton  principle  for  irreversible  processes  in  two- 
component  media 


4.4.1  The  governing  equations 


Let  us  consider  irreversible  processes  in  two  component  media.  In  order  to  include  into  con¬ 
sideration  the  Darcy  and  Navier-Stokes  dissipation  we  assume  the  presence  of  the  viscous  mass 
force  Pa  vie  viscous  stress  tensor  Also  we  extend  the  Hamilton  principle  and 

formulate  it  in  the  following  nonholonomic  form 


^  r  dt  I  dopL-^  r  dtj  doY,  (pl 

dr  Jto  Jn  Jto  Jn  a  ^ 


DgUu 
'  Dt 


—  V  -Pg^i'a 

‘‘^a.vis  ^  3 


)  = 


(4.25) 
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Copying  the  derivation  of  the  system  (4.24)  we  arrive  at  the  following  generalization  of  the 
governing  system 

~  Pb^i^  ~  P  "I"  Lpb^iPb  + 

+Vi{pphLp^)  -  "^jipLvjukb^'kib)  +  H.vis  +  (4-26) 

This  reduces  to  Eq.  (4.24)  for  conservative  systems. 

4.4.2  Energy  relations  and  general  properties 

For  each  constituent,  a,  define  the  elastic  stress  tensor 

”  ^P^"^ jUka^kia  +  PPaLp^Sf  ,  (4.27) 

and  the  “interaction”  force,  or  momentum  supply  in  the  terminology  of  mixture  theory  [41,  42], 

Ga  —  Pa^il'  ~  P  ^Pa^iPa  »  (4.28) 

then  the  equilibrium  conditions  (4.26)  become 

Wjf  =  Vi  {Ft  +  Pt^,)  +  Gi  +  PL,  ■  (4-29) 

Note  that  the  sum  of  the  interaction  forces  at  a  point  is  zero, 

E^i  =  0-  (4-30) 

a 

The  sum  of  the  elastic  and  viscous  stresses,  Pt-\-Pa.m»^  can  be  considered  the  total  stress  in  phase 
a.  We  caution  that  these  deiinitions  are  not  unique  because  the  vectors  G\  could  be  changed  by 
adding  divergence-like  terms. 

Define  the  total  energy  density  per  unit  mass  of  the  system, 

^  =  («i) 

and  define  the  flux  vector  for  component  a 

n  =  -<pL  .  (4-32) 

where 

=  (4-33) 

Pb 

is  the  mean  velocity.  Then  it  follows  from  Eq.  (4.26),  with  the  details  in  the  Appendix,  that 
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This  is  the  pointwise  energy  dissipation  relation.  It  is  interesting  to  note  that  the  flux  vector  of 
phase  a  contains  a  term  proportional  to  the  relative  velocity  of  a  with  respect  to  the  mean  flow. 

Integrating  the  bcdance  law  (4.34)  over  the  domain  Q  and  using  Eq.  (4.15)  leads  us  to  the 
following  identity 

5  i,  ^  =  /n 


which  is  the  energy-dissipation  relation  for  the  volume  Cl.  In  the  absence  of  dissipation  we  can 
establish  the  energy  conservation  equation 


dt 


=  0. 


(4.36) 


Alternatively,  let  C  be  the  Lagrangian  density  per  luiit  volume, 

C=z  pL, 


(4.37) 


in  terms  of  which  the  governing  equations  (4.26)  become 


Ph^i^pb  ^ji^VjUbb-^kib) 

0 


+PL,  +  ^iPiL-  (4.38) 

This  form  of  the  equilibrium  equations  is  perhaps  simpler.  It  is  interesting  to  note  that  if  C  is 
additive  in  the  following  sense, 

^  ^  Pa  I^ai^Pa'i  ^ j'^ka)  >  (4.39) 

a 

then  the  equations  (4.38)  reduce  to 

{pbPb,pbS"^  -  PbLbyjutb^'kih)  +  Pb.vi»  +  ^jPLi,-  (4.40) 

These  equations  are  coupled  only  through  the  viscous  terms.  Thus,  an  additive  volumetric 
Lagrangian  density  leads  to  decoupling. 

4.4.3  Dissipation  effects:  Constitutive  relations 

For  the  stresses  we  Umit  ourselves  with  the  traditionail  choice  of  the  rheology 

fL.  =  (4.41) 

where  the  tensor  of  viscous  coefficients  is  symmetric  with  respect  to  the  first  and  second 
pcdrs  of  indices  and  to  interchcinge  of  these  pairs.  Also,  based  on  Eq.  (4.34),  it  should  satisfy  the 
foUowmg  dissipation  inequality 


(4.42) 
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In  order  to  establish  a  constitutive  relation  for  the  viscous  body  force  P^^^g  we  impose  the 
following  “natureil”  demands: 

a)  dissipativity 

b)  additivity 

=  (4.44) 

a 

c)  the  “action-reaction”  equality 

FL=-fL,  (4.45) 


d)  linearity  in  the  velocities 


abc^J*^  ' 


(4.46) 


The  constraint  a)  is  a  natural  consequence  of  the  energy  identity  (4.34),  while  the  conditions  b) 
and  c)  imply  that 

=  (4-47) 

b 

which  is  similar  to  the  general  result  (4.30). 

We  now  examine  the  implications  of  these  constr2unts  for  two-component  media.  First,  Eq. 
(4.45)  implies  that 


i^i  =  =  0, 

Combining  (4.46)  and  (4.48)  we  get 


jpt  _  jpt 

^12  —  “-^21  • 


4"  7?222^j2  —  0> 


(4.48) 


(^121  +  -^211)  +  (^122  +  -^212)  —  0. 

Equations  (4.49)  lead  us  to  the  following  formulais  of  the  coefficients 

Fni  —  -^112  =  -^221  =  -^222  — )  0 


(4.49) 


•^121  —  ~-^211  -  Xl  )  .£^122  -  “-^212  — 


(4.50) 
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Combining  Eqs.  (4.48)  -  (4.50)  we  obtadn  the  relationships 

^12  “  “-^21  ~  ^122^i2  =  Xl'^jl  +  X2'^j2i 

ui  ^  jpi  T>i  ^  jpi  jpi  _  pi 

^l^vis  -^21 J  ■^2.vi8  “^12  -^21  ’ 

Hence, 

Plvis^il  +  Plvi»Vi2  =  -  f’i2)  =  -{XlVjx  +  X'ivj2){Vii  -  Vis)  ■ 

The  inequality  (4.43)  will  be  guaranteed  provided  x-l  =  —Xs  =  a  positively  definite  matrix. 
Indeed,  in  this  case  we  get 

+  Plvi,^i2  =  -  Vi2)iVjl  -  Vjs)  <  0. 

With  the  above  mentioned  choice  we  obtain  the  following  constitutive  equation: 

PLi.  =  -D'^ivii  -  Vis),  PLu  =  -D'\vi2  -  Vii),  PLis  =  -D'^ivia  -  Vib).  (4.51) 

Inserting  the  constitutive  relations  (4,41)  and  (4.51)  in  (4.26)  we  arrive  at  the  following  master 
system: 


Dt  ph^  jUkb^  +  ^iiPPbLpf^) 

-^jipLvju,,Akii>)  -  D'^ivib  -  vi^)  +  V,-  (z?£VfcU,5)  ,  (4.52) 

or  equivalently, 

i^Vl^kVib)  ,  (4.53) 

where  the  interaction  force  is 

^6  “  “^a  ~  P^  I^pa^iPa  + 

-Pa  +  Ppb^iPb  +  LVjUkb^i^jV’kb^  ■  (4.54) 


4.5  Examples 

4.5.1  One-component  models. 


Example  1.  Nonlinear  elastic  solid. 
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The  governing  system  of  a  nonlinear  one-component  solid  can  be  derived  from  the  master  system 
(4.53)  if  we  neglect  viscous  forces  and  choose  the  Lagrangian 


L  =  T  -  E  =  -v^vi  -  E{VjUk). 

In  this  case  P^^  of  (4.27)  is  the  Cauchy  stress  tensor 

=  p^^Ai , 
dVjUu 


(4.55) 


and  there  are  no  viscous  or  interaction  terms.  The  equilibrium  equations  (4.53)  are  simply 

(4-66) 


Example  2.  The  Navier- Stokes  fluid. 

The  governing  system  of  a  nonlineair  one- component  viscous  fluid  can  be  derived  from  the  master 
system  (4.52)  if  we  ignore  the  drag  force  and  choose  the  Lagrangian  and  the  viscosity  tensor 
in  the  following  form: 

L^T-E  =  iu’vi  -  E{p),  +  v,  ,  (4.57) 

where  are  the  volume  eind  shear  viscosities,  respectively.  The  elastic  stress  is  hydrostatic, 

defined  by  a  pressure  p. 


nii  cH  2^^ 

P^'  =  -p6^\  P  =  P-^, 

(4.58) 

and  the  viscous  stresses  P’i.  Eire 

vta 

# 

+  I/.  +  V4’)  . 

(4.59) 

The  equilibrium  equations  (4.53)  are  now 

# 

Dv^ 

p—  =  -V‘p  + v,pv  . 

(4.60) 

As  cdways,  when  combined  with  the  mass  conservation  equation  (4.3)  the  system  (4.59),  (4.60) 
is  closed  in  terms  of  p,  v*. 

4.5.2  Two-component  models. 

Let  us  assume  that  the  kinetic  energy  T  is  a  quadratic  form  of  the  velocities 


(4.61) 
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Using  (4.61)  we  get 

0  =  5  E  =  j  E  "i- 

b 

In  the  case  of  isotropic  two-phase  medium  we  get 

iT'i  =  icabS^^,  D'^  =  AS'^. 

Combining  (4,62)  and  (4.63)  we  get 

^  =  ^kacVla, 

where  the  symmetric  matrix  kac  is  defined  as 

kac  ~  “h  kcc^  . 


(4.62) 


(4.63) 


(4.64) 


Example  3.  Liquid  two-component  mixture. 

We  obtain  the  simplest  possible  model  of  a  two-component  liquid  mixture  by  choosing  the  La- 
grangian  and  viscosity  tensors  in  the  following  form: 

L  =  T  —  E  ^  V  ^  kfiijV^Vii)  “  E(pa), 

ab 


The  partial  stresses  are  hydrostatic, 


Pb=-Pb6^\  Pb  =  PPb 


dE 

dpb' 


(4.65) 


(4.66) 


and  the  interaction  forces  are  given  by 


^b  —  Pb^  i  ^2  ^ac^aVic  kqfyVjg^  jV^  pEp^^  ip}) 

—  '^ic  {j^b^ca^ a  Pakcb^ "b  Pa^pf,^ iPb  PbEp^^iPa  • 


(4.67) 


Inserting  these  into  the  master  equations  (4.53)  we  obtain  the  following  system  of  equilibrium 
equations  for  the  two-liquid  mixture 


Pb 


Dt 


I  ^  ^  k(i})Via  I  —  ^  ^  '^ic  (^Pbkca^ Pakcb^ i'^})^  H“  Pa^pi^ ipb  Pb^Pa^iPo, 

\P^  a  /  c 

-ViippbE,,)  -  Aivib  -  Via)  +  V,-  [u^bSiVkvt  +  u,b  {Vivi  +  V^vij,)]  .  (4.68) 
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Example  4«  Solid  two-component  mixture. 


The  simplest  model  of  a  two- component  solid  mixture  corresponds  to  a  Lagrangian  of  the  form: 

L  =  T  -  E  =  -  ^  kabvlvib  -  E{ViUja).  (4.69) 

^  ab 

Inserting  (4.69)  in  (4.53)  we  obteiin  the  following  master  system  for  a  two-component  solid  elastic 
mixture: 

Pb^  -Pb^^i  +  pY,^<^bVjb^ivl  =  -pb^iE 

+pEv,u,,Vi^jUkb  +  ^i{pEviu,,Akib)  -  D'^[vib  -  Vi,)  +  V,-  {DfJ,,VkVia)  .  (4.70) 


Example  5.  Fluid  permeable  elastic  solid. 

We  choose  subscript  s  for  the  solid  and  /  for  the  liquid  component.  The  simplest  model  of 
two-component  solid  mixture  is  obtcuned  by  choosing  the  Lagrangian  in  the  following  form: 

L  =  T  ^  E  =  "1“  (4.71) 

Inserting  (4.71)  in  (4.53)  yields  the  following  master  system: 

P»^  -f  A!,/v,/))  =  p,  {kf,Vjg  +  kffVjf)  VjvJ  -  Pf  {kttVj,  +  k^fVjf) 

-p,ViE  +  pE^jui^yiVjUk,  +  ^  jipE^jut.A'kis) 

-D'^Vi.-vjf)  +  Vj{Di^!:lVkVu),  (4.72) 

Pf^  (*//^*/  +  ^»/^'«)  j  =  Pf  -f  ktfVjf)  Vivl  -  p,  (Aj/.Vj,  -I-  kffVjf)  Viv^ 

~Pf^iF  +  pEp^^ipf  —  ^i[ppfEp^) 

-  Vi,)  +  Vj  {Dy':!,,VkVif)  .  (4.73) 


4.6  Conclusion 

We  have  established  a  gener^d  framework  for  dealing  with  fully  nonlinear  problems  for  two- 
component  media.  The  central  result  is  Eq.  (4.26)  which  gives  the  equilibrium  equations  in 
Eulerian  coordinates  for  an  arbitrary  Lagrangian  density.  These  governing  equations  are  a  direct 
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consequence  of  the  generalized  Hamilton’s  principle  in  the  presence  of  dissipative  effects,  f.e.,  the 
Onsager-Sedov  principle  of  Eq.  (4.25).  We  have  illustrated  the  application  of  the  general  govern¬ 
ing  equations  to  specific  examples  of  single  and  two-phase  continua.  In  particular,  the  standard 
theories  of  solids  and  fluids  are  included.  The  full  potentieil  of  the  theory  is  its  applicability  to 
multi-phase  continua. 

4.7  Appendix:  Some  details 

Here  we  derive  the  energy  identity  (4.34).  We  start  with 

b 


EDa^  Da  I k  T 


(^)  -  .  (A.1) 


We  next  simplify  each  of  these  terms  in  turn.  First,  using  (4.17), 

Da  (  dL  \  J,  Df)  dL  r-v  k  (  l  J  ^  r7 

\pbDh  (  p  dL\  dL  pbDb  p  (  I  I 


= 

b 


p  Dt  \pbdvt )  dvt  p  Dtpb 


I  i  i\  t-i 


which  implies  that 


u  Da  I  dL  \ 

a  b  \  o  /  b 


Mr. 

bPh 


Di  (  p  dL\  dL  Dh  p 


Dt  \  pb  dv^  I  dv^  Dt  pi 


dL 


(A-2) 


and  is  defined  in  (4.33).  By  making  use  of  the  governing  system  (4.26)  we  get  the  following 
result: 

-^3  +  Vi{ppbLp^)  +  Pligjb  +  ■  (•^•2) 


Use  of  Eq.  (4.3)  gives  us 

Db  P 
'Dtpb 


Db  P  DbPa  P  DbPb  r  .  T7  ..I  (..l  ..l\T7  .  I 

Pb-fT,-  =  l^  —  --—  =  2^[-pa^l^a-iv.-Vb)'^lPa  + 


Dt  Pb  Dt 
=  V,(p(t;[-F'))  , 
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implying 

=  (A.4) 

Combining  Eqs.  (A. 2),  (A.3),  and  (A.4)  gives  us  the  first  sum  in  the  left  member  of  Eq.  (A.l), 

-Vj  (pLvjun,Akii,^  +  ^iipPbLpJ  +  P^igh  +  '^jPyit,b 

+V/(pi„.(F'-t;'))]  .  (A.5) 

Now  consider  the  other  terms  in  (A.l).  Using  the  mass  conservation  Eq.  (4.3)  again  gives 

a  b  a  b 

=  (A.6) 

b 

Also,  by  using  of  the  last  of  the  Eq,  (4.18)  (with  the  change  of  the  parameter  r  for  t)  we  get  the 
following  equation: 

=  HYj  [Aj^b  -  (v^  -  u^)  Vfc  Vj«/b] 

a  b  ^  b 

=  ^pLvjuit  [AirfcVjVb  -  Ub  VfcVjUib  +  U*'VfcVjUjb]  .  (A.7) 

b 

By  combining  (A.l),  (A.5),  (A.6),  and  (A.7)  we  obtain  the  relation  (4.34),  where  the  flux  vectors 
are  defined  in  (4.32). 


Chapter  5 

Acoustoelasticity  of  fluid-saturated 
porous  media 


The  general  theory  for  small  dynamic  motion  superimposed  upon  large  static  deformation,  or 
acoustoelasticity,  is  developed  for  isotropic  fluid-filled  poroelastic  solids.  Formulae  are  obtained 
for  the  change  in  acoustic  wave  speeds  for  arbitrary  loading,  both  on  the  frame  and  the  pore 
fluid.  Specific  experiments  are  proposed  to  find  the  complete  set  of  third  order  elastic  moduli  for 
an  isotropic  poroelastic  medium.  Because  of  the  larger  number  of  third  order  moduli  involved,  7 
as  compared  with  3  for  a  simple  elastic  medium,  experiments  combining  open-pore,  closed-pore, 
jacketed,  and  unjacketed  configurations  are  required.  The  details  for  each  type  of  loading  are 
presented,  and  a  set  of  possible  experiments  is  discussed.  The  present  theory  is  applicable  to 
fluid-saturated,  bi-connected  porous  solids,  such  as  sandstones  or  consolidated  gremular  media. 


5.1  Introduction 

The  hnear  theory  of  poroelastic  fluid  saturated  media  is  a  mature  subject,  with  the  classic  studies 
of  Frenkel  [6]  and  especizdly  of  Biot  [7]  standing  as  the  signiflc2int  achievements  in  its  modem 
development.  Several  different  methods  now  exist  for  arriving  at  governing  equations,  including 
the  theory  of  interpenetrating  continua  [8],  and  the  method  of  “homogenization”  based  on  two- 
sczile  asymptotic  expainsions  [59,  10].  Despite  the  fact  that  these  approaches  do  not  always  yield 
the  seime  governing  equations,  eind  further  general  theoretic^  studies  are  required,  they  have 
certainly  proved  useful  in  severed  branches  of  civil  etnd  geo-engineering  and  the  Earth  sciences  [11]- 
[13],  and  also  in  unexpected  fields,  such  as  low  temperature  physics  [14]  and  pattern  formation 
in  polymer  gels  [15].  However,  there  etre  definite  limitations  to  the  linear  theory.  For  example, 
a  proper  ainalysis  of  the  dynamics  of  large  zunplitude  sound  in  sediments,  or  of  small  amplitude 
sound  in  rocks  tmder  laurge  confining  stress  eu:e  obviously  impossible  within  the  framework  of  a 
linpar  theory.  Even  seismic  waves  of  very  small  intensity  at  a  large  distance  from  their  epicenter 
should  be  studied  in  the  framework  of  a  nonhneeir  theory  because  of  the  weU-known  effect  of 
the  accumulation  of  nonlinear  distortions  leading  to  the  “gradient  cateistrophe”  phenomenon 
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There  has  been  renewed  interest  in  the  nonlinear  acoustics  of  rocks,  based  partly  on  several 
observations  of  distinct  nonlinear  effects.  These  include,  the  change  of  the  in  situ  velocities  of 
seismic  waves  [18]“[20],  and  direct  measurement  of  harmonic  distortion  [21,  22].  These  exper¬ 
iments  have  clearly  demonstrated  that  the  in  situ  nonlinear  elastic  moduli  of  rocks,  soil  and 
sediments  are  much  greater  than  the  linear  ones.  The  nonlinear  behavior  of  fluid-filled  poroe- 
lastic  solids  is  not  as  well  understood,  and  further  progress  depends  upon  accurate  experimental 
measurement  of  the  effective  nonlinear  moduli  of  poroelastic  media. 


The  determination  of  the  linear  moduli  has  been  the  subject  of  numerous  experimental  and 
theoretical  studies  [23,  24,  12].  An  isotropic  poroelastic  medium  has  four  static  moduli,  which 
may  be  measured  by  a  combination  of  stress-strain  experiments  on  “jacketed”  and  “unjacketed” 
samples.  The  role  of  the  jacketing  is  to  constrain  either  the  fluid  pressure  or  its  mass.  The 
nonlinear  moduli,  in  contrast,  should  be  determined  by  measurement  of  essenticJly  nonlinear 
effects.  One  of  the  simplest  nonlineetr  phenomena  is  the  elastoacoustic  effect,  whereby  the  speeds 
of  smeiU  Eimplitude  waves  are  changed  by  applying  stress  or  stredn.  The  theory  of  acoustoelasticity 
has  been  thoroughly  discussed  for  purely  elastic  materials  [25]- [27],  and  is  now  commonly  used 
in  ultrasonics  [28]- [30],  and  has  been  adapted  to  multiphase  materials  [5]. 


The  subject  of  this  Chapter  is  the  acoustoelastic  effect  in  the  context  of  poroelasticity.  The 
theory  is  based  on  a  nonlinear  generalization  of  the  classic  Biot  theory.  Biot  himself  discussed 
thermodynamic  aspects  of  a  nonlinear  theory  of  poroelastic  media,  but  did  not  specify  any 
nonlinear  stress-strain  behavior  (see  Papers  15,  16,  18  and  19  in  ref.  [7]).  There  has  been  some 
work  on  nonlinear  poroelasticity  within  the  framework  of  the  theory  of  interpenetrating  continua 
[8],  and  cdso  using  the  method  of  two-sceJe  homogenization  [31,  2].  Three  sources  of  nonlineetrity 
are  traditionally  distinguished  in  elasticity:  (i)  the  physical  nonlinearity,  that  is  the  nonlinearity 
of  the  constitutive  relations  relating  the  stress  and  the  displacement  gradient;  (ii)  nonlinearity 
of  the  universal  equations,  such  as  the  equations  of  conservation  of  mass,  momentum,  etc.;  (iii) 
geometric  nonlinearity,  which  results  from  the  nonlinear  relationship  between  the  deformation 
gradient  and  the  tensor  of  finite  deformations. 


In  order  to  keep  the  analysis  as  simple  as  possible,  in  this  Chapter  we  ignore  aU  dissipative 
effects  and  concentrate  on  the  first  kind  of  nonlinearity  only  because  it  appears  to  be  the  most 
significant  for  geophysical  materials.  We  begin  with  a  discussion  of  acceleration  waves  in  nonlinear 
poroelastic  media.  These  are  exact  solutions  independent  of  the  state  of  prestress  in  the  medium, 
and  they  reduce  to  the  weU-known  Biot  fast  etnd  slow  waves  for  unstressed  isotropic  materials. 
These  general  results  serve  as  the  basis  for  considering  the  principal  body  waves  propagating 
in  a  slightly  prestressed  isotropic  nonlinear  saturated  poroelastic  substance.  Finally,  we  discuss 
application  of  the  general  formulae  for  the  stress  dependence  of  the  wave  speeds  to  severed 
specifically  designed  experimental  configurations. 
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5.2  Nonlinear  poroelasticity 

5.2*1  Governing  equations 

All  further  analysis  is  based  on  the  following  governing  equations  of  the  nonviscous,  fluid-fiUed 
poroelastic  medium: 


where 


^  dt^  dt^  ~ 

(A.la) 

dt^  dt^  ~  ^  P’3' 

(A.lb) 

iJ  \  ^ 

“  duij  dC 

(A.2) 

Here  u*  is  the  displacement  of  the  solid  skeleton,  w*  is  the  relative  displacement  of  the  fluid, 
are  the  stresses,  p  is  the  fluid  pressure,  W’(C,Um,n)  is  the  poro-elastic  potential,  p  and  pf  are 
the  averaged  density  and  the  density  of  the  fluid,  respectively;  —(  =  is  a  divergence  of  the 
relative  displacement  of  the  fluid,  is  the  “instantaneous”  magnitude  of  a  symmetric  tensor  of 
permeability  (the  value  of  the  permeability  hereditary  operator  at  t  =  0).  Also,  are  the  spatial 
coordinates,  the  Latin  indices  take  the  vedues  1,  2,  and  3;  summation  over  repeated  indices  is 
implied,  and  a  comma  followed  by  a  Latin  suffix  symbolizes  partial  differentiation.  We  refer  to 
the  publications  [7],  [9],  [31],  and  pcirticularly  [2]  and  for  the  motivation  behind  Eqs.  (A.l)  and 
(A.2). 

By  inserting  the  vedue  of  (  in  terms  of  w*  we  can  rewrite  the  equations  of  motion  as 


d^w' 


fd^uj 


(A.3a) 


(A.3b) 


In  these  and  subsequent  equations  we  use  the  following  notation  for  the  derivatives  of  the 
potential  energy  function  W:  =  dW/duij,  =  d^WjdC^  =  d^W/duijdC,  = 

d^W/duijduk,i,  etc. 


5.2.2  Acceleration  waves 

An  acceleration  wave-front  is  a  propagating  surface  defined  such  that  displacements  ^lnd  their 
first  derivatives  are  continuous  across  this  surface,  but  the  second  and  higher  derivatives  cuid,  in 
peirticuleir,  the  accelerations  possess  finite  jumps.  Acceleration  waves  are  the  most  convenient 
object  of  theoretical  study  for  both  linear  and  nonhnezir  dynamics.  They  have  been  investigated 
thoroughly  by  Hadamard  [32]  and  Thomas  [33],  while  Chen  [34]  provides  a  more  recent  review  for 
elastic  materials.  Continuity  of  the  first  derivatives  imposes  geometric  and  kinematic  constraints 
on  possible  jumps  of  the  second  derivatives  since  the  “tangential”  components  of  the  second 
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derivatives  are  continuous.  These  constraints  are  known  as  compatibility  conditions,  etnd  are  es- 
sentiaiUy  geometric  conditions,  originally  developed  by  outstanding  geometers  like  Hadamard  [32], 
Levi-Civita  [35],  and  Thomas  [33]  (fortunately  all  three  were  outstcinding  mathematical  physi¬ 
cists  and  their  works  are  accessible  to  physicists).  In  particulcir,  the  jumps  of  second  derivatives 
across  acceleration  wave-fronts  satisfy  [34,  36] 

+ 

=  hiC^y  (A.4a) 

-H 

=  Hic\  (A.4b) 

where  hk  =  [ui^pqltvPn*^  and  Hk  =  are  the  second  order  amplitude  vectors  of 

discontinuity;  is  the  unit  normal  to  the  wave-front,  and  c  is  the  velocity  of  the  front.  The 
compatibility  conditions  indicate  that  the  two  vector  functions  hk  and  Hk  completely  define  the 
jumps  of  60  partial  derivatives  of  the  displacements. 

In  order  to  determine  the  acceleration  vectors  themselves  one  has  to  extract  some  additional 
dynamic  information  from  the  governing  Eqs.  (A.3).  These  equations  are  not  valid  at  the  wave- 
front  since  the  second  derivatives  are  undefined  at  the  front.  By  definition,  only  one-sided 
second  derivatives  are  defined  at  the  front,  and  these  one-sided  limits  satisfy  Eqs.  (A.3).  We  first 
subtract  term-wise  Eq.  (A.3a)  for  the  two  one  sided  limits  teiking  into  accoimt  continuity  of  the 
first  derivatives  ctnd,  then  do  the  same  operation  with  Eq.  (A. 3b).  Then,  evaluating  the  jumps 
of  the  second  derivatives  using  the  second  order  compatibility  conditions  of  Eq.  (A.4)  we  get  the 
following  linear  cdgebraic  system  for  the  amplitude  vectors 

(^pc'^6^'‘-W'^'^^njni)hk+(pfC^S^'‘  +  Wi^njn'‘)Hk  =  0,  (A.5a) 

[pfcH^'^  +  W^‘nin^)hk+{pfC^Kit-Wan'n'‘)Hk  =  0,  (A.5b) 

where  is  the  inverse  of  the  permeability  matrix  ,  which  is  assumed  to  be  invertible. 
The  pressure  p  is  continuous  across  the  wave  front,  but  its  gradient  is  not,  specifically 

[p,i]l  =  qrii 

At  the  same  time,  it  foUows  from  the  Eqs.  (A.2)  and  (A.4)  that 

[p,i]l  =  W^'ninihk  -  WccTiin’^Hk.  (A.7) 

Combining  the  previous  two  equations  allows  us  to  eliminate  the  quantity  ri^Hk  in  favor  of  the 
jump  in  the  pressure  gradient,  g.  Then  Eq.  (A. 5b)  provides  an  expression  for  ZT*, 


(A.6) 


=  HiUjUk, 


d^Wi 


dt^ 


""  hiTljTlky 


d^ui 


dt^ 
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After  some  manipulation,  Eqs.  (A.5)  imply  that 


where 


-W^^^^njni)hk-  K^^rijq  =  0, 

(A.9a) 

T^kl  L  ( K'^niUj  1  \ 

K  nihk  -j-  I  2  ~  TIT-  )  9  — 

(A.9b) 

=  pS*'^  -  PfK^'^, 

(A.lOa) 

(A.lOb) 

j^ik  ^  /r**' 4- wf 

(A.lOc) 

Equations  (A. 9)  form  a  system  of  4  linear  cdgebraic  equations  with  respect  to  the  4  unknowns 
hk  and  q.  This  system  has  a  non-zero  solution  when  the  determinant  vanishes,  leading  to  a  4-th 
order  polynomial  equation  in  c^.  The  coefficients  of  the  polynomial  depend  on  Uk,i,  C  and  the 
orientation  of  the  unit  normal  n*,,  and  thus,  the  same  is  true  for  the  velocities  c.  Each  real  root 
generates  ein  associated  non-zero  solution  (hjt.g).  Further  discussion  of  the  4-th  order  system 
can  be  found  in  [37].  On  the  other  hjoid,  the  system  defined  by  Eq.  (A.5)  is  6-th  order,  because 
it  involves  the  unknowns  hi  and  Hj.  The  4-th  and  6-th  order  systems  are  cormected  by  the  fact 
that  the  latter  possesses  the  root  =  0  of  multiplicity  2.  Thus,  when  c^  =  0  the  system  of  Eq. 
(A.5)  has  infinitely  many  solutions  of  the  form  hk  =  0,  Hk  =  Rk  where  Rk  is  an  arbitrary  vector 
in  the  2D-space  orthogonal  to  the  wave  normed  nt .  In  summmy,  the  propagating  waves  aind  the 
non-zero  wave  speeds  C2in  be  foimd  from  either  the  4-th  order  or  the  6-th  order  systems. 


5.2.3  Principal  directions 

A  principal  direction  is  defined  as  one  in  which  an  eigenvector,  such  as  h*,  is  aligned  with  or 
orthogonal  to  the  wave  normal  n*.  For  example,  consider  a  longitudinal  wave  in  the  direction  of 
ni,  that  is,  a  solution  of  the  form 

hk  —  ,  Hk  —  HjXk’  (^'^^) 

Substituting  from  Eq.  (A.ll)  into  Eq.  (A.5)  and  contracting  with  nj  gives 

[pc^--W'^^^ninjnkni)h-t[pfC^  +  Wi‘nini)H  =  0,  (A.12a) 

(^PfC^-\-Wi‘nini)h-\-[pfC^Kil^,nink--Wc()H  =  0.  (A.12b) 

It  follows  immediately  from  Eqs.  (A. 12)  that  the  corresponding  eigenvalues  are  defined  by 
a  biquadratic  cheiracteristic  equation,  whereas  the  amplitudes  h  and  H  are  connected  by  the 
relation 

pc^  —  W^^^^niTijUkTii  _ 

Pf<^^K^Lnink  -  l^CC 


(A.13) 
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A  transverse  wave,  on  the  other  hand,  satisfies  hk  =  hck,  Hk  =  Hck  where  e*  is  a  polarization 
vector  with  n*'ej.  =  0  and  =  1.  Under  these  circumstances,  Eqs.  (A.5)  become 

ekh pfC^e^H  =  0,  (A.14a) 

[pfcH^'^-[-W^%v^)ekh-^PfC^KilekH  =  0,  (A.14b) 

implying  that  the  transverse  wave  velocity  is  given  by 

C?  =  {p-pf  {Ki^v^peg)~^)  ^  W'^'^^eiekTijUi .  (A.15) 

When  the  material  is  isotropic  then  all  directions  Eire  principal  directions,  i.e.,  all  acceleration 
waves  Eire  either  longitudinEil  or  trsuisverse  in  nature. 

5.3  The  acoustoelastic  effect 

5.3.1  Third  order  moduli  for  poroelasticity 

We  now  consider  the  effects  of  prestress  on  the  longitudinal  and  transverse  waves  in  sm  isotropic 
medium.  The  permeability  for  an  isotropic  poroelastic  medium  is  =  k  while  the  potentiEd 
W  is  a  function  of  (  Euid  the  principal  invariants  of  the  symmetric  strain  Sij  =  {ui 
Thus, 


W  =  W{Iuh,h,0,  (A.16) 

where 

h  =  l2  =  l  ejgtkr.  h  =  (A.17) 

and  is  the  third  order  Edternating  tensor.  To  within  third  order  terms  in  Sij  Eind  (  the 
potential  W  can  be  approximated  by  the  polynomial 

=  rn/i'  +  r2/2  +  ric/iC  +  r(cC'  +  rin/f  +  ri2/i/2 

+r3/3  +  Ticc/iC'  +  Tiic/f  c  +  r2c/2C  +  rc«c"-  (a.is) 


We  consider  the  following  field  of  perturbations  with  respect  to  the  reference  (unstressed) 
configuration 


6ui  =  I  a^njX^Tii  +  ^  (A.19) 

V  >1=1.2  / 

where  eu  Emd  e2i  are  two  mutusdly  orthogonEil  transverse  vectors,  i.e.,  both  Eu:e  orthogonal  to 
the  n*  direction.  The  perturbation  magnitude  is  defined  by  e,  with  |£|  <C  1  by  assumption,  and 
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the  parzaneters  a^f,  ua,  A  =  1,2,  and  A  are  0(1)  aind  independent.  The  prestrain  is,  from  Eq. 
(A.19), 

eij=  jQ7^njni+  ^  e.  (A.20) 

\  A=1,2  J 

Differentiating  Eq.  (A. 18)  in  the  vicinity  of  the  reference  configuration  (e  =  0)  we  obtain  the 
second  order  elasticity  tensors  in  terms  of  the  elastic  moduli: 

+  /X  ,  wi^  =  -aM  6^^,  W^c  =  M,  (A.21) 

where  Ac,  a,  and  M  are  the  Biot  parameters  [38].  In  fact,  Biot  defines  the  linearized  stress-strain 
relations  in  terms  of  a  quadratic  energy  potential  similar  to  Eq.  (A. 18).  Comparison  of  the 
quadratic  terms  in  the  latter  with  equation  (3.4)  of  Biot  [38]  immediately  implies  that 

Ac  =  2rii-fr2,  p  =  ®  =  “riij/(2r(^^).  (A.22) 

There  axe  meiny  different  notations  used  for  describing  the  lineeir  response  of  saturated  porous 
media,  and  even  Biot’s  own  notation  evolved  from  the  time  of  his  early  work  on  consohdation 
(see  Paper  1  of  ref.  [7])  to  his  later  work.  His  1962  paper  [38]  provides  a  good  comparison  of  the 
notations  used. 

For  the  third  order  quantities  we  get 

=  (Pac  +  2rn(:)  i  Fac  , 

Wii  =  Wccc  =  GTiCO  (A.23) 

^n^klmn  ^  ^  SijSklSmn  +  ^  Pj  -f  +  /Irn^kn  ^ 

_iri2  -I-  6^^  -1-  6""^  . 

In  the  absence  of  poroelastic  effects,  the  third  order  elastic  moduli  Wijkimn  are  equivalent  to  the 
standard  third  order  moduli  Cijkimm  or  Cjjk  in  the  concise  Voigt  notation  [28],  In  order  to 
compare  the  moduli  here  with  the  more  common  notation,  we  note  the  equivalences 

Till  =  ^12  =  2  (^112  C'lii) ,  “  2  ^  3(7ii2  +  2C123)  ♦  (A. 24) 

The  connection  with  other  notations  for  third  order  moduli  can  be  inferred  from  the  table  in 
Green’s  review  [28]  which  compares  many  different  systems  of  notation.  We  note  for  future 
reference  the  following  formulae 

W^^^^^^rtinjUkTii  =  6TniS^^  +  2ri2  -  n^n^) ,  (A.25a) 

Pxa^”*”  -  ^Pa  (^’""  -  n’”n’‘ -  e”*e”)  .  (A.25b) 

iU  4U 
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5.3.2  Wave  speed  dependence  on  the  applied  strain 

For  each  deformed  configuration  there  are  three  orthogon^ll  directions  of  principal  straiin,  zmd  the 
smcdl  eimplitude  waves  which  propagate  in  these  directions  su'e  either  longitudinal  or  transverse. 
We  will  now  derive  approximative  formulae  for  the  incremental  cheinges  of  the  velocities  of  these 
principal  waves  to  leading  order  in  e. 

First,  the  veiriation  in  the  velocity  of  the  transverse  wave  with  the  polarization  vector  ea 
follows  from  Eqs.  (A.15)  and  (A.20)  as 

=  (fi- 

-  (p- (w'^^‘'”"eciecknjniem„  +  W^^'°‘eciecknjnieA^  .  (A.26) 

Note  that  the  densities  p  and  pf  are  unchanged  because  we  are  using  the  reference  or  Lagrangian 
description,  as  opposed  to  the  current  or  Eiderian  description.  Conservation  of  mass  requires 
that  these  densities  are  constant.  The  structure  constant  k  is  also  assumed  to  remain  the  same 
imder  the  deformation.  Combining  Eqs.  (A.20),  (A. 25)  and  (A.26)  we  arrive  at  the  formula  for 
the  elasto- acoustic  effect  for  trcuis verse  waves, 

“  2(p  -  Kpf) 

In  peurticular,  the  split  in  the  velocities  of  two  transverse  waves  with  the  same  propagation 
direction  but  different  poleirizations  is,  from  Eq.  (A.27), 

Scj,  -  eVs 

aA  -  OLB  2(p  -  Kpf) ' 

In  order  to  estabhsh  analogous  formulae  for  longitudinal  principal  waves  we  first  write  the 
system  of  Eq.  (A. 5)  as 


(A.28) 


((ai  +  02  +  (xpj)  Ti2  +  (oi  +  02  —  ac)  Ta  +  T2C^)  •  (A.27) 


AX  =  CX , 


(A.29) 


where  X  =  (h,  and  the  symmetric  matrices  A,  C  are 


—W^niUi 

P 

PS 

A  = 

,  C  = 

-W^UiUi 

.  Pf 

! 

1 _ 

(A.30) 


The  inertial  matrix  C  is  again  unchanged  by  the  initicd  deformation,  and  the  new  wave  speeds  de¬ 
pend  upon  the  variation  in  A.  By  mziking  use  of  Eqs.  (A. 18)  through  (A. 23)  we  can  approximate 
this  to  within  first  order  in  £, 


A  =  A^  +  £AS 


(A.31) 
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where 


Ac  +  2p 

aM 

6Pni 

-2Piic  ■ 

A"  = 

aM 

M 

,  =  a^ 

-2Pii^ 

2Picc 

(A.32) 


6Piii  +  2Pi2 

-2Pii^  -  P2f 

2P11C 

-2P1CC 

+  A 

— 2Piic  -  P2C 

2Pi« 

.  -2ricc 

srccc 

+  (Ofl  +  Ot2) 


The  eigenvalue  problem  for  the  matrix  A°  was  first  studied  by  Biot  in  1956  (see  Paper  8  of 
ref.  [7])  who  established  the  existence  of  two  different  longitudinal  modes:  the  “fast”  and  “slow” 
waves  with  velocities  cj  and  Cg,  respectively.  The  eigenvectors  Xy  and  X,  associated  with  these 
velocities  can  be  calculated  by  making  use  of  the  following  equation  implied  by  Eq.  (A. 13), 


(A.33) 


—  M 

Corrections  to  the  Biot  fast  and  slow  velocities  due  to  the  incremental  deformations  C2in  now  be 
calculated  using  standard  perturbation  theory,  leading  to  the  formula: 

AhX^X^^ 


Sc^  =  £ 


_ 


(A.34) 


Equations  (A.27),  (A.32),  and  (A.34)  show  that  by  measurement  of  elasto-acoustic  effects  we  can, 
in  principle,  find  the  third  order  elastic  moduli  by  solving  a  system  of  linear  algebraic  equations. 
In  the  next  section  we  will  examine  some  possible  experimental  configurations  with  this  purpose 
in  mind. 


5.4  The  elasto-acoustic  effect  in  experiments 

5.4.1  Experimental  nomenclature  and  moduli 

The  apphcation  of  the  previous  ideas  to  a  poroelastic  sample  requires  the  ability  to  impose  various 

types  of  stress  states  ifedl  seven  third  order  moduli  are  to  be  determined,  i.e.,  {Pin,  ri2,  Ps,  Pif^,  Piio  r2C> 

We  first  discuss  the  types  of  measurements  employed  in  determining  the  four  second  order  moduli 

{Pill  P2>  PiC>  Tcc)  or  {-^c)  o,  M}.  The  second  order  moduli  eire  related  to  one  2mother  via 

Eq.  (A.22),  and  yield  the  classic  Biot  theory  for  linear  poroelasticity,  for  which  [38] 

Tjj  =  AgU  6ij  -|-  -j-  uy  j)  otM^  6ij  , 


p  =  -aMu^k-^-M^. 


(A.35a) 

(A.35b) 


Note  that  p  is  the  pore  pressure,  and  Tij  is  the  total  or  confining  stress.  The  effective  stress, 
Tij  +  ap%,  is  also  commoidy  used.  For  example,  if  a  hydrostatic  confining  stress  =  -Pc^ij 
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is  applied,  then  the  effective  pressure  is  Pe  =  Pc  -  otp^  The  main  point  is  that  there  are  two 
independent  stress  variables,  or  pressures  if  the  system  is  hydrostatically  loaded.  Table  1  hsts 
and  defines  three  distinct  bulk  moduli  (inverse  of  compressibility),  Kc,  and  Km^  associated 
with  hydrostatic  deformation  of  undrained,  jacketed,  oxid  unjacketed  Scimples,  respectively.  We 
also  define  the  bidk  modulus  for  the  fluid  alone,  Kf^  eind  the  “coefficient  of  fluid  content”  7 
[23,  38], 


=a(1 _ l-\ 

Kf  Vf  dp’  ^  Km}' 

in  terms  of  which  the  Gassmeinn  relation  [38]  is, 

_  1  a 


(A.36) 


(A.37) 


Different  types  of  experiments  have  been  used  to  determine  the  second  order  elasticities 
of  fltud-saturated  poroelastic  media.  We  shall  discuss  them  briefly  in  order  to  avoid  possible 
confusion  in  terminology.  The  terms  “open”  and  “closed”,  “jacketed”  and  “unjacketed”,  and 
“drained”  and  “undrained”  are  all  relevant  to  the  problem,  but  have  different  connotations.  We 
wiU  next  describe  the  four  states  of  deformation  illustrated  in  Figure  5.1,  and  then  consider  their 
mathematical  description,  and  their  implications  for  elas to- acoustic  measurements. 

(a)  The  “closed-pore”  system  corresponds  to  constcincy  of  the  fluid  content  or  mass,  and 
in  the  geometrically  hnear  approximation  this  implies  that  ^  =  0.  The  closed-pore  experiment 
is  the  simplest  for  theoretical  amalysis  but  is  rather  difficult  to  realize  in  practice,  at  least  for 
a  high  level  of  wetting.  Conceptually,  it  requires  that  the  sample  is  covered  by  an  impervious 
closed  deformable  “jacket”,  see  Figure  5.1a.  (b)  Alternatively,  the  “open”  porous  system,  in  the 
terminology  of  Biot  and  Wilhs  [23],  corresponds  to  the  case  of  constant  fluid  pressure,  or  p  =  0. 
However,  it  is  more  consistent  to  refer  to  this  as  a  “drzdned”  test  and  to  treat  it  as  a  special 
case  of  a  more  general  “open”  configuration  with  p  =  const.  Such  an  experiment  is  sketched 
in  Figure  5.1b  and  is  regarded  as  the  “jacketed”  test.  The  amount  of  absorbed  liquid  is  not 
fixed  in  this  experiment  since  it  can  leave  the  sample  through  the  tube  connected  with  the  ah. 
The  “unjacketed”  tests  are  sketched  in  Figure  5.1c  and  Figure  5. Id.  Although  the  system  as  a 
whole  is  again  contciined  within  an  impervious  jacket,  the  porous  solid  sample  can  now  exchange 
fluid  with  the  surroimding  free  fluid  in  the  reservoir,  (c)  A  “conventional”  “imjacketed”  test  is 
shown  in  Figure  5.1c.  In  this  case  the  hydrostatic  stress  on  the  solid  is  caused  by  pressurizing  the 
surrounding  fltiid,  and  therefore  both  the  skeleton  and  the  fluid  are  hydrostaticsdly  stressed  with 
the  same  pressure,  Tij  =  —Pc^ij  and  p  =  pc-  This  type  of  conventional  “unjacketed”  test  does 
not  permit  tri-axial  deformation  of  the  sample,  and  is  thus  too  restrictive  for  the  acoustoelastic 
effect,  (d)  The  experiment  sketched  in  Figure  5. Id  might  be  helpful  in  avoiding  this  drawback, 
although  it  is  presumably  rather  difficult  to  realize  practically.  In  what  follows  we  call  it  a  “tri- 
axicil  imjacketed”  test,  and  it  is  cheiracterized  by  a  constant  pressure  p  in  the  pore  fluid,  and  a 
tri-axial  state  of  stress  in  the  skeleton,  implying  three  independent  elements  for  r^j. 

The  four  states  outlined  above  are  as  follows.  First,  for  the  “closed-pore”  jacketed  test  we 
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have  (see  Table  1  for  definitions) 


0 

II 

^  p^-Pi 

K,  K' 

(A.38) 

In  the  “open-pore”  jacketed  test,  on  the 

other  hand. 

0 

II 

II 

1 

p 

(A.39) 

In  the  “conventional”  nnjacketed  test  the  skeleton  is  in  a  state  of  hydrostatic  stress,  experiencing 
the  same  pressure  as  the  fluid,  ie.^pc  =  p.  In  the  case  at  hand  the  skeleton  deformation  is  purely 
dilatational,  and  Eqs.  (A. 35)  then  imply 

Suij  =  -  Sij  ,  SC  =  TP,  (A.40) 

where  7  is  defined  in  Eq.  (A. 36).  Finally,  in  the  tri-axial  “unjacketed”  test  both  the  fluid  pressure 
and  the  deformation  of  the  skeleton  are  simultaneously  kept  under  control.  In  this  case  we  can 
express  the  fluid  dilatation  parameter  in  terms  of  these  as 

K  =  (A.41) 


5.4.2  Wave  speeds  for  specific  experimental  configurations 


The  above  formulae  together  imply  specific  relations  for  the  change  in  the  velocities  of  the  three 
wave  types:  transverse,  and  fast  ctnd  slow  longitudinal  waves.  For  the  first  three  cases  below,  the 
loading  is  hydrostatic,  defined  by  the  confining  pressure  pc.  The  change  in  the  transverse  speed 
is  therefore  independent  of  its  polarization  direction. 

(a)  The  “closed-pore  jacketed”  tests:  for  the  transverse  wave, 

dc^  _  ^12  +  iFa 

dpc  “  2Kc{p-  KpfY 

and  for  the  fast  eind  slow  waves, 

dc2 

dpc  -  CijXo^Xo^  ’ 


where 

2  — fiFiii  — |ri2  2rii^ -h  |r2(; 

[  2rnc  +  ir2c  -2T,cc 

(b)  The  “open-pore  jacketed”  tests: 

dcf  _  ri2  -f  |r3  -f  ar2^ 
dpc~  2K{p-Kpf) 


(A.44) 


(A.45) 
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The  variation  in  the  fast  and  slow  wave  speeds  is  given  by  Eq.  (A. 43),  where  now 

j  “SFiii  ”  3  ri2  -  2arii^  +  I  r2C  +  2ari^^ 

2rii^  +  |r2^  +  2ari^^ 

(c)  The  conventional  “unjacketed”  tests: 

dc^  _  ri2  +  \Tz  -  7-^Mr2C 
dpc  2Km{p  -  /cp/) 

The  fast  and  slow  speed  chcinges  are  given  by  Eq.  (A.43)  with 


(A.46) 


(A.47) 


B 


1 

Km 


— 6riii  —  |ri2  -  ^iKm^uc 


2Tii(  +  ^T2(  +  2jKm^i(( 


2riic  +  3  r2c  +  27iirMricc 


.  (A.48) 


(d)  The  “unjacketed”  tri-cixial  tests:  The  three  principal  strains  aj\r,  ai,  a25  and  the  pore  pressure 
p  =  eP  are  now  independent  parameters,  and,  according  to  Eq.  (A.41), 


A  —  a  {aj^  +  ai  +  02)  +  — . 

The  change  in  the  three  wave  speeds  are  then  given  by  the  general  formulae  in  Eqs. 
(A.32)  and  (A.34). 


(A.49) 

(A.27), 


5.5  Discussion  and  conclusions 

The  acoustoelastic  effect  is  the  simplest  nonlinear  elastic  effect  which  can  be  used  to  determine 
nonlinectr  elastic  moduli.  To  this  end  one  should  measure  changes  of  velocities  of  acoustic  waves 
induced  by  prestrain.  Although  it  is  of  higher  order  as  compared  to  linear  effects,  acoustoelasticity 
can  be  reliably  detected  in  sedimentary  materials  since  the  nonlinear  elastic  moduli  are  of  several 
orders  higher  that  their  linear  moduli  [21].  The  experiments  can  be  performed  either  on  core 
samples  or  in  situ;  for  instance,  by  experimental  measurement  of  strain  induced  changes  in 
seismic  waves  velocities  due  to  the  solid  Earth  tides  [18].  In  fact,  it  is  found  that  the  in  situ 
nonlinearity  is  as  large  or  larger  them  observed  in  the  laboratory,  and  may  be  due  to  the  Ijtrge 
scale  heterogeneity  in  the  Earth  [18].  Such  measurements  can  provide  valuable  information  about 
the  Ecurth’s  interior,  just  as  measurements  of  the  velocities  themselves  allow  one  to  determine 
linecir  moduli. 

In  the  present  Chapter  we  have  established  relations  for  the  strain  induced  changes  of  acoustic 
wave  speeds  in  fluid-filled  poroelastic  media.  We  have  concentrated  on  the  longitudinal  and 
trains  verse  waves  in  isotropic  poroelastic  substances,  which  become  slightly  anisotropic  under 
triaxial  prestradn.  These  changes  depend  not  only  on  the  prestrain  of  the  solid  skeleton  but  also 
on  the  drainage  regime  of  the  fluid  component.  Thus,  the  acoustoelastic  effect  is  different  for 
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drained,  undrained,  jacketed,  and  unjacketed  tests,  and  explicit  formulae  for  each  configuration 
are  described. 

The  7  nonlinear  elastic  moduli  characterizing  an  isotropic  poroelastic  substance  can  be  found, 
in  principle,  by  a  sequence  of  acoustoelastic  tests  and  using  the  formulae  presented  here.  Thus, 
the  nonlinecu-  elastic  modulus  Ta  may  be  determined  directly  using  the  rather  simple  formula  of 
Eq.  (A.28)  for  the  split  in  the  velocities  of  two  transverse  waves  having  the  same  direction  of 
propagation  but  different  polarizations.  Then,  using  the  formulae  in  Eqs.  (A.42)  and  (A.45)  one 
cm  determine  the  nonlinear  moduli  ri2  md  12^  by  making  shear  wave  measurements  imder  uni¬ 
form  (hydrostatic)  prestrain  with  “closed-pore  jacketed”  and  “open-pore  imjacketed”  conditions. 
The  same  information  cm  be  extracted  using  shear  wave  data  from  a  conventional  “imjacketed” 
test  md  Eq.  (A.47).  Compatibility  of  the  numerical  results  for  the  moduli  determined  through 
different  experiments  cm  be  used  as  the  indication  of  acceptability  of  the  model  of  m  isotropic 
poroelastic  substmce.  The  determination  of  the  remaining  4  nonlinear  moduli  cm  be  based  on 
measurements  of  the  “fast”  md  “slow”  longitudinal  waves,  using  Eqs.  (A.34),  (A.43),  (A.44), 
(A.46)  md  (A.48).  Again,  the  fact  that  there  are  more  relations  thm  moduli  allows  one  to  check 
the  admissibility  of  the  model  of  m  isotropic  poroelastic  medium. 

I^ble  I 


Compressibility 
-  definition 

In  current 
notation 

Terminology,  notation  £ind  references 

1  _  -1  av 

Kc  -  V  dpc 

C 

1 

“closed”  compressibility,  [38] 
undrained  compressibility,  c„  [39] 

1 

open  or  jacketed  compressibility,  [38,  23] 
drained  or  matrix  compressibility,  c  [39] 

1  -iav\ 

Km  ~  V  dp^  \j 

»=Pc 

m 

gredn  or  pore  compressibility,  [40],  Cg  [39] 
imjacketed  compressibility,  6  [23] 

Table  1.  Different  compressibilities  of  a  porous  medium,  with  definitions  of  Kc,  K  md  Km- 
Here,  p  md  pc  are  the  pore  md  confining  pressure,  respectively.  Note  that  Ac,  p,  a  md  M  are 
defined  in  Eq.  (A. 22),  md  A  =  Ac  —  a^M  [38].  The  relations  between  these  md  other  moduli 
are  explored  in  greater  detail  by  Kiimpel  [39],  md  edso  by  Brown  md  Korringa  [40]. 


Chapter  6 


Waves  and  flutter  in  fluid  conveying 
elastic  media 


A  fluid  flows  through  a  periodically  layered  system  comprising  alternating  elastic  solid  and  liquid 
constituents.  An  exact  analysis  is  performed  to  consider  the  existence  and  stability  of  small 
acoustic  waves  and  disturbances.  The  presence  of  the  flow  introduces  the  possibility  of  flow 
induced  flutter.  Unstable  waves  are  generadly  possible  for  M  of  order  unity,  where  M  is  the  shear 
wave  Mach  number,  but  appear  for  much  lower  veJues  for  asymmetric  flexural  type  motion.  In 
that  case  it  is  foimd  that  a  critical  wavenumber  exists,  indicating  that  the  system  is  inherently 
unstable  to  long  wavelength  disturbances. 

6.1  Introduction 

Fluid-saturated  permeable  elastic  sohds  exhibit  distinct  phenomena  not  seen  in  either  the  sohd 
or  the  liquid  phase.  The  most  interesting  acoustic  feature  is  the  appearance  of  the  Biot  “slow” 
wave,  predicted  by  Biot  (1956)  and  observed  by  Fiona  (1980)  in  water  saturated  sintered  beads. 
Similar  acoustic  effects  have  since  been  seen  in  other  systems,  including  layered  periodic  solid- 
hquid  configurations  (Fiona  et  al.  1987).  The  periodicity  of  the  layered  structure  permits  a 
more  precise  theoretical  analysis  than  for  a  disordered  system,  eind  very  good  agreement  has 
been  observed  between  theory  (Rytov  1956,  Schoenberg  1983)  and  experiment  (Fiona  et  al. 
1987).  Acoustics  of  such  layered  solid-liquid  structures  is  of  great  interest  in  various  applications 
relating  to  physics,  mechanical  engineering,  Earth  science,  etc.  (Brekhovskikh  1981),  because  the 
structure  displays  the  essential  features  of  a  realistic  permeable  elastic  medium.  In  all  previous 
studies  of  this  system  it  was  assumed  that  both  the  solid  cmd  fluid  constituents  are  at  rest  in  the 
“ground”  equilibrium  configuration.  On  the  other  hand,  it  is  weU-known  that  sufficiently  fast 
fluid  motion  through  a  deformable  solid  structure  Cctn  destabilize  the  system  causing  flutter-like 
phenomena  (Bolotin  1963).  So  fzir  flutter  has  been  studied  using  engineering  theories  of  plates 
and  shells  for  the  solid  phase,  combined  with  hydrodynamic  theories  for  the  liquid,  which  are 
sometimes  based  on  the  hypothesis  of  plane  cross-sections,  rather  than  with  the  help  of  much 
more  rehable  theories  of  bulk  hydrodynamics  and  elasticity  (Fung  1955,  Bolotin  1963). 
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The  purpose  of  this  Chapter  is  to  report  on  new  results  2ind  physiccil  phenomena  caused  by 
the  non-smcill  relative  velocity  of  the  fluid  propagating  between  solid  elastic  layers.  Some  of 
these  effects  are  a  direct  consequence  of  the  flow,  and  vanish  in  the  equilibrium,  no-flow  state. 
The  study  of  acoustic  distiirbEinces  and  flutter-like  phenomena  in  the  presence  of  flow  requires 
a  nonlinear  basis,  and  hence,  our  study  begins  with  an  exact  nonlinear  formulation.  We  then 
linearize  the  equations  and  interface  boundary  conditions  in  the  vicinity  of  a  stationary  state  in 
which  parallel  isotropic  elastic  layers  are  undeformed,  while  the  fluid  moves  with  the  constant 
velocity  V^,  Our  objective  is  to  explore  the  dependence  on  of  the  dispersion  equations  of 
acoustic  waves  in  the  small  vicinity  of  the  steady,  uniform  flow  configuration. 

6.2  Formulation  and  general  solution 

6.2.1  Nonlinear  equations 

Let  us  consider  a  nonviscous  layered  media  consisting  of  alternating  solid  and  fluid  layers.  We 
choose  the  Euleri^ln  description  of  continua  in  order  to  simplify  exact  formulation  of  the  interface 
boundary  conditions.  We  emphasize  that  the  relative  displacements  of  the  fluid  and  solid  con¬ 
stituents  are  not  small,  generzJly  speaking.  We  denote  by  **  the  Eulerian  Cartesian  coordinates, 
and  assume  that  the  coordinate  is  directed  along  the  undisturbed  elastic  layer.  The  bulk 
equations  within  the  solid  and  fluid  constituents  and  the  interface  boundary  conditions  eire  then 


+  (A.la) 

v:  =  ^  +  vivjui,  (A.ib) 

/dvi  .\ 

=  -Vp,.  (A.1C) 

^  +  V,(wl7)  =  0,  (A.ld) 

priij  =  (A.l«) 

V/AT,  =  VfNj,  (A.lf) 


where  ps  and  pf  aie  the  actual  densities  of  the  constituents,  V*  and  V}  are  their  velocities,  Ul 
is  the  displacement  of  the  sohd,  P/*  and  pf  are  the  Cauchy  stress  tensor  of  the  solid  and  the 
fluid  pressure,  respectively,  and  is  the  unit  normeil  of  the  interface  (which,  generedly  speaking, 
changes). 

The  system  of  equations  (A.l)  permits  a  stationeiry  solution  in  which  the  sohd  layers  are  at 
rest  and  undeformed  whereas  the  fluid  moves  with  constcuit  velocity  V^  parallel  to  the  elastic 
layers.  Linearizing  the  system  of  equations  (A.l)  we  arrive  at  the  governing  system  for  the  small 
disturbzinces. 
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6.2.2  Linearization 

We  consider  small  dynamic  motion  superimposed  on  the  state  of  uniform  pressure  and  flow  in 
the  fluid,  corresponding  to  fluid  velocity  V°,  fluid  pressure  p°,  tind  densities  and  The  solid 

is  also  assumed  to  be  in  a  state  of  imiform  static  initial  stress,  P/'  =  Pi*,  either  hydrostatic 
(p/*  =  -p°6^*^  or  otherwise,  with  the  initial  deformation  homogeneous  and  defined  by  U*  =  P*. 
Thus,  we  let 


V} 

=  W'  +  uj, 

(A.2a) 

Pi 

=  P°  +  Py 

(A.2b) 

ui 

= 

(A.2c) 

PV 

(A.2d) 

where  and  cdl  the  remaining  dyn2tmic  quantities  («’ ,  V,*,  Pt—pl,Pf—Pf,  and  p)  aie  small. 
We  consider  motion  in  the  plane,  where  x^  is  the  coordinate  in  the  layering  direction. 

Then  equations  (A.l)  implies  the  following  linearized  equations 


(A.3a) 

II 

(A.3b) 

n  /  dv)  dv\  \ 

(A.3c) 

(A.3d) 

I 

II 

b 

(A.3e) 

(A.3f) 

The  equations  for  the  solid  phase,  (A.3a,b),  must  be  supplemented  by 
relation  for  the  small  stress  of  the  form 

a  lineeir  stress-strain 

(Tji  =  Cjifc/  V 

(A.4) 

Therefore,  eqs.  (A. 3a)  and  (A.3b)  imply  that  P*  satisfies  the  usual  equations  of  linear  dynamic 
elasticity. 
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The  fluid  equations  (A.3c,d)  require  2in  additional  linear  equation  of  state, 

Pf-  P)  = 

where  Cf  is  the  speed  of  sound.  When  combined  with  eqs.  (A. 3c)  and  (A. 3d),  this  implies  a 
convective  wave  equation  for  the  small  pressure, 


6.2.3  Traveling  wave  solutions 

Let  us  consider  “traveling- wave” -solutions  of  the  lineeurized  system,  i.e.,  the  solutions  proportional 
to  There  are  two  sorts  of  such  solutions  which  we  caU  symmetric  and  asymmetric, 

respectively.  For  the  symmetric  mode  the  shape  of  each  layer  remains  symmetric  with  respect 
to  its  median.  The  dispersion  equation  for  this  mode  is  derived  below  as 


/  tanh(fc6g.)  ..  .  I  P)  (.  ^c^anhCfe^g.) 

y  cl)  c  ) 


(A.8) 


Here  we  use  the  following  notation:  c  —  wikis  the  velocity  of  the  traveling  wave,  2Hg  and  2Hf 
are  the  equilibrium  thicknesses  of  the  fluid  and  solid  layers,  respectively,  =  (l  - 

=  (l  —  and  =  ^1  -  (c-  ^  ,  where  Ct  and  c/  are  the  velocities  of  bulk 

transverse  and  longitudinal  waves,  respectively,  within  the  undeformed  solid  layer,  and  cy  is  the 
bulk  sound  velocity  within  the  undisturbed  fluid.  Also,  and  pj  are  the  imdisturbed  densities. 

In  the  asymmetric  mode  opposite  edge  points  of  each  layer  have  the  same  verticcd  velocity. 
The  dispersion  equation  for  the  asymmetric  mode  is  the  following: 


/  tanh(fe6g,) 

\  cl)  ^  pt\  c  )  if  ct  tanh(A:6if,) 


(A.9) 


6.2.4  Derivation 


The  dispersion  relations  can  be  derived  using  impedance- type  concepts.  Consider  the  layer  of 
solid  of  thickness  subject  to  a  normal  stress  on  either  face  such  that 

the  motion  is  either  symmetric  or  asymmetric.  The  shear  stress  on  both  faces  is  zero.  Similarly, 
consider  the  moving  fluid  layer  subject  to  a  pressure  disturbance  p  =  po  ^  again  either 

symmetric  or  asymmetric.  The  symmetry  or  asymmetry  means  we  need  only  consider  half  of 
the  unit  period  of  the  system,  that  is,  the  solid  and  fluid  half  layers  in  <  0  and 

0  <  <  ify,  respectively.  For  each  half-layer  problem,  let  Vg  be  the  vertical  velocity 

Vy  or  at  the  stuface  =  0.  Define  the  effective  impedcUices, 


a)  _  ^0  „(ol)  _  PO 


(A.10) 
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where  a  =  ±1  indicates  the  symmetry.  Thus,  a  =  1  eind  —1  correspond  to  the  symmetric  and 
asymmetric  configurations,  respectively. 

In  the  absence  of  the  flow  the  dispersion  equations  for  guided  waves  of  the  system  are  defined 
by  the  identity  =  0,  which  follows  from  the  continuity  of  traction  and  normed  velocity 


at  the  interface.  That  is,  from  eqs.  (A.3e)  eind  (A.3f)  with  V°  =  0, 

=  -p  and  V/  =  Vy  at  =  0. 


(A.11) 


However,  in  the  presence  of  flow  the  kinematic  constraint  of  equation  (A.ll)2  must  he  replaced 
by  (A.3f).  For  the  pcirticular  dependence  this  simplifies  to 


vt=  1 


V.  at  =  0. 


Thus,  the  dispersion  relation  in  the  presence  of  flow  becomes 


+  (l  -  ^)  =  0,  a  =  ±  1. 


A  simple  calculation  based  on  equations  (A. 7)  and  (A. 3c)  shows  that 
^(±1)  =  Hf)fK 


(A.12) 


(A.13) 


(A.14) 


For  simplicity,  we  assume  the  solid  layer  to  be  isotropic,  in  which  case  standard  analysis  gives 


=  1 4^  ( UiEf  (cotUklfH.))^^  -  (2 


mt  icoth{k^tHs)f^  -  2  -  -2  {coth{k^iH,)) 


(A.15) 


6.3  Asymptotic  limits 

6.3.1  Short  wavelength  asymptotics:  two  half  spaces 

We  now  consider  several  asymptotic  cases  of  the  general  residts  defined  by  equations  (A. 8)  and 
(A.9).  In  the  short  wavelength  limit  \kHg\  1*^/1  >>  1  the  dispersion  equations  for  both 
modes  lead  to  the  same  equation: 


(-8’— 


(A.16) 


By  introducing  nondimensional  p2u:8uneters  M  /ct,  s  =  c/ct,  si  =  cifct,  Sf  =  Cf/ct,  and 
S  =  we  CEUi  rewrite  the  last  equation  in  the  following  form: 
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This  is  identiceil  with  the  Rayleigh  equation  for  Pf  =  0  (Rayleigh  1885,  Achenbach  1973).  ^lnd 

it  gives  the  Scholte  equation  for  V°  =  0  (Scholte  1948,  1949).  By  solving  it  numerically  one 

can  find  that  the  roots  for  the  interface  wave  speed  c  =  u^/k  become  complex  and  the  system  ^ 

becomes  unstable  for  V°/ct  =0(1).  For  example.  Figure  6.1  shows  the  merging  of  two  real  roots 

as  M  is  increased  from  zero.  In  this  example,  complex  roots  appejir  for  V"  greater  than  about 

1.793  Ct. 


6.3.2  Long  wavelength  asymptotics:  symmetric  mode 

The  long  wavelength  asymptotics  \kHg\  ~  \kHf\  <<  1  appears  to  be  quite  different  for  the 
syimnetric  and  asymmetric  modes.  For  the  sjrmmetric  mode  we  get 


or  in  non-dimensional  form 


(2_.2)^_4(i-^  +^s2(.-M)^ 


/  1-4  \ 


1 

A 


=  0, 


(A.18) 


(A.19) 


where  h  =  Hf/Hg,  This  can  be  rewritten  as  a  fourth  order  polynomial  equation  for  5, 

((.  -  Mf  -  ,})  (.’  -  .5)  ^  +  (,  -  Mf  -  .?)  =  0,  (A.20) 

where  4>  =  -I-  Hg)  is  the  porosity,  and  Sp  =  Cpjct  where  Cp  is  the  speed  of  a  longitudinal 

“plate”  wave;  thus,  =  2/(1  -  v).  We  note  that  the  symmetric  waves  for  the  long  wavelength 
asymptotic  limit  of  an  isolated  plate  in  a  fluid,  \kHg\  «  1,  \kHf\  oo,  are  contained  in  equation 
(A.20)  as  the  limiting  case  >  1. 

At  M  =  0  we  get  a  quadratic  equation  with  respect  to  which  has  two  physiceilly  meaningful 
roots:  the  velocities  of  the  so-called  Biot  “fast”  and  “slow”  waves  (Rytov  1956,  Schoenberg 
1983,  1984,  Fiona  et  al.  1987)  The  velocities  of  both  waves  do  not  depend  on  the  direction 
of  propagation.  The  non-zero  relative  velocity  V°  of  the  fluid  destroys  the  equivalence  of  the 
opposite  directions.  For  small  M  the  last  equation  allows  one  to  And  the  magnitude  of  the 
velocity  splitting  for  each  type  of  wave.  Thus,  let  sq  >  0  he  the  non-dimensional  wave  speed 
relative  to  Ct  for  M  =  0,  fast  or  slow,  eind  let  si  be  the  other  speed,  slow  or  fast,  then  the  split 
velocities  are  ±so  +  ^M,  where 


The  Biot  fast  and  slow  wave  roots  for  M  =  0  simplify  when  both  constituents  Me  incom¬ 
pressible,  that  is  si,Sf  00,  Sp  =  2.  The  fast  speed  then  becomes  infinite  but  the  slow  roots 
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are 


c 


V^S 
6  +  h^ 


(A.22) 


6.3.3  Asymmetric  modes  -  flexural  waves 


In  the  asymptotic  limit  of  \kHs\  <<  1  the  dispersion  equation  (A. 9)  for  the  asymmetric  mode 
gives 


s 


2 


=  0. 


(A.23) 


Flexural  type  waves  are  found  using  the  scaling  c  =  Ct  kHg  r,  where  r  is  a  parameter  of  order 
unity.  Substituting  this  into  equation  (A.23)  yields, 


r 


2 


(’• 


M  y  tanh(fcg/e/)  0 

kHj  kHfCf 


(A.24) 


The  flexural  wave  dispersion  relation  for  a  plate  in  vacuo  follows  from  this  with  ^  =  0  as  c  =  ±ct 
where  Cf,  =  Ct  \kHs\  v^2/3(l  -  u)  is  the  phase  speed  of  a  flexural  or  bending  wave  on  a  plate. 
The  dispersion  relation  for  a  fluid-loaded  plate  in  the  absence  of  flow  is  obtained  from  (A.24)  by 
setting  M  =  0.  We  have  retained  the  parameter  (f  in  (A.24)  rather  than  set  it  to  unity  in  order 
to  be  consistent  with  steindard  analyses  for  fluid  loaded  plates,  e.g.  Junger  and  Feit  (1986). 

If  the  fluid  layer  is  also  very  thin  compared  with  the  wavelength,  i.e.  \kHf\  «  1,  then 
equation  (A.24)  reduces  to  a  quadratic  equation  in  r,  yielding 


where 


My/3  _  W 

Cf) 


(A.26) 


When  M  =  0,  the  roots  yield  c  =  ic^/x/l +Th,  which  correspond  to  flexural  waves  on  cin  isolated 
plate  which  has  the  bending  stiffness  of  a  single  elastic  layer,  and  the  mass  of  a  single  period  of 
the  solid-liquid  system.  That  is,  the  effect  of  the  fluid  is  just  an  added  mass  as  it  moves  in  phase 
with  the  flexuTcd  motion.  For  small  values  of  M,  or  equivalently  e,  the  two  roots  are  regular 
perturbations  of  the  flexural  wave  roots.  For  large  €,  on  the  other  hand,  we  have  the  possibility 
of  of  two  complex  roots  when  the  discriminant  of  equation  (A. 25)  goes  to  zero.  .The  existence  of 
complex-conjugate  roots  indicates  that  the  flow  causes  a  flutter-like  instability.  This  occurs  for 
€  >  €c  where  the  critical  value  is 


ee  =  . 


(A.27) 
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In  summary,  anti- symmetric  disturbances  of  the  layered  system  are  imstable  for  very  long 
(e  >  1  or  \kHs\  <C  M)  wavelength,  but  short  wavelengths  (e  <C  1)  are  stable.  There  is  a  critical 
wavelength  above  which  ciU  disturbances  are  unstable,  and  it  is  defined  by  e  =  Cc,  or  A:  = 
where 


1 

2 


(A.28) 


Finally,  in  order  to  analyze  flutter  of  an  isolated  elastic  layer  we  let  \kHf\  oo  in  the 
dispersion  equation  (A. 23)  for  the  asymmetric  mode.  Using  the  same  variables  as  before,  we 
consequently  obtain  the  following  equation: 


,  HkH.f 

3(1 -k)'^ 

When  the  fluid  is  incompressible  this  equation  becomes  a  quadratic  with  roots 

V^S 


c  = 


8-^\kH.\  6  +  \kH,\ 


11  + 


^(1  - 
\kH.\ 


incompressible  fluid. 


(A.29) 


(A.30) 


The  possibility  of  complex  conjugate  roots  again  shows  that  flutter  instability  occurs  for  long 
wavelength  perturbation  of  the  system.  That  is,  the  system  is  stable  (unstable)  for  |Aj|  >  kc 
(1*1  <  kc)  where  the  flnite  wavenumber  which  defines  the  onset  of  the  flutter  regime  is 


and  A  is  the  unique  positive  root  of 


A"  +  A2  = 


3^02 


(A.31) 


(A.32) 


The  case  of  an  isolated  elastic  plate  in  a  mean  flow  has  received  extensive  coverage  in  the  literature 
using  thin  plate  equations  (Kirchhoff  theory).  For  example,  Crighton  eind  Oswell  (1991)  discussed 
the  response  of  a  thin  plate  in  an  incompressible  flow  to  a  line  drive.  They  iilso  provide  a  useful 
review  of  the  literature  on  the  stability  of  flexible  surfaces  and  compliemt  coatings  in  the  presence 
of  mean  flow. 


6.4  Summary 

The  flow  of  compressible  fluid  through  a  layered  medium  provides  a  very  rich  system  for  studying 
the  phenomenon  of  acoustics  in  fluid  conveying  structures.  Starting  from  the  exact  nonlinear 
equations  of  motion  we  have  derived  the  equations  for  small  dynamic  disturbances  superimposed 
upon  a  steady  flow  configuration.  The  possible  wave  types  for  the  periodically  layered  medium 
may  be  distinguished  by  their  pairity,  symmetric  and  asymmetric,  each  of  which  displays  quite 
distinct  stability  characteristics. 
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The  general  dispersion  relation  for  symmetric  modes  is  given  by  equation  (A. 8),  with  short 
cind  long  wavelength  limits  in  equations  (A. 16)  and  (A. 18),  respectively.  It  is  found  that  insta¬ 
bility  is  possible  only  for  flow  speed  on  the  order  of  the  bulk  wave  speeds,  i.e.,  for  M  =0(1). 

Regarding  asymmetric  modes,  the  genercil  dispersion  relation  of  equation  (A. 9)  has  the  same 
short  wavelength  limit  as  that  for  symmetric  modes,  equation  (A. 16).  Thus,  disturbances  of  short 
wavelength  become  imstable  only  for  M  =0(1).  However,  disturbances  of  very  long  wavelength 
are  potentially  unstable  in  the  presence  of  flow  for  both  the  periodic  system,  as  indicated  by 
equation  (A. 25),  and  for  an  isolated  plate,  from  equation  (A. 30).  The  associated  long  wavelength 
wave  types  are  analogous  to  flexureJ  waves  on  plates  in  vacuo.  For  a  given  flow  speed,  or  M, 
there  is  a  critical  wavenumber,  key  such  that  quasi-flexural  waves  are  stable  for  k  >  key  but 
instability  is  possible  for  all  k  <  ke. 


Chapter  7 
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Figure  6.1  The  real  roots  for  flow  over  a  solid  half  space,  from  equation  (A.17).  For  simplicity, 
we  have  taken  the  fluid  emd  solid  as  incompressible  {si,Sf  oo),  and  ^  =  1.  The  emergence  of 
complex  roots  occurs  at  M  «  1.793  for  this  case. 
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INTRODUCTION 

Biot’s  theory^’^  for  linear  dynamics  of  fluid-filled  porous 
media  is  by  now  well  established,  with  several  experimental 
investigations  that  corroborate  it  very  well.  For  example,  the 
data  of  Fiona  et  al.^  display  both  the  fast  and  slow  wave  in  a 
layered,  alternating  fluid/solid  system.  It  has  recently  been 
emphasized  that  sandstone  and  other  porous  granular 
media*^’^  display  strongly  nonlinear  behavior,  as  compared 
with,  for  example,  water  or  metals.  The  nonlinearity  param¬ 
eter  B/A  for  water,  which  is  the  ratio  of  a  third-order  elastic 
modulus  to  the  second-order  bulk  modulus,  is  approximately 
5.  The  analogous  quantity  for  sandstone  can  be  on  the  order 
of  10'^.'^  What  can  one  expect  for  the  nonlinear  behavior  of  a 
strongly  nonlinear  porous  sandstone  saturated  with  water? 
The  nonlinear  acoustical  properties  of  the  constituents  are 
disparate,  and  it  is  not  clear  how  this  will  affect  the  nonlinear 
acoustics  of  the  fast  and  slow  waves  of  the  poroelastic  me¬ 
dium.  This  paper  is  a  first  step  toward  an  understanding  of 
the  nonlinear  mechanics  and  dynamics  of  porous  media,  par¬ 
ticularly  the  interplay  between  the  nonlinearities  of  the  con¬ 
stituents.  We  focus  on  the  problem  of  deriving  the  governing 
equations,  the  strain  energy  functions,  and  describe  some  lin¬ 
ear-wave  solutions  for  a  layered  medium. 

Our  objective  is  the  governing  nonlinear  equations  for  a 
poroelastic  medium — the  nonlinear  generalization  of  the 
Biot  model.^  We  employ  the  “two-scale”  technique  of  ho¬ 
mogenization  for  heterogeneous  media  with  disparate  length 
scales,  which  has  been  used  to  obtain  the  linear  Biot 
theory. Homogenization  leads  directly  to  the  macroscopic 
equations,  and  defines  the  microproblems  which  uniquely 
determine  the  coefficients  in  them.  Burridge  and  Kostek^^ 
recently  derived  the  nonlinear  poroelasticity  equations  for  a 
system  with  a  granular  solid  skeleton  whose  elastic  response 
is  governed  primarily  by  deformation  at  grain  contacts.  This 
is  assumed  to  dominate  the  nonlinear  effects,  to  the  extent 
that  the  fluid  can  b^  considered  as  linearly  elastic.  Unlike 
Burridge  and  Kostek,^^  we  do  not  assume  that  the  solid  non¬ 
linearity  overwhelms  the  fluid’s,  but  keep  both  on  an  equal 


footing  with  the  purpose  of  comparing  their  interaction. 

Three  sources  of  nonlinearity  are  traditionally  distin¬ 
guished  in  continuum  mechanics.  First,  there  is  the  physical 
nonlinearity,  that  is,  the  nonlinearity  associated  with  the  de¬ 
pendence  of  the  stress  tensor  on  the  tensor  of  finite  deforma¬ 
tions.  The  second  one  is  the  nonlinearity  of  the  universal 
equations  of  mass,  momentum,  etc.,  that  is,  those  equations 
which  have  validity  for  all  specific  models  of  the  continuous 
medium.  The  third  one  is  geometrical  nonlinearity,  for  ex¬ 
ample,  the  nonlinear  relationship  between  displacements  and 
the  tensor  of  finite  deformation.  This  classification  of  nonlin¬ 
earities  is  not  absolute,  and  some  changes  occur  when  one 
switches  from,  for  instance,  the  Eulerian  to  the  Lagrangian 
description  (e.g.,  the  universal  momentum  equation  appears 
to  be  nonlinear  in  the  Eulerian  description  and  linear  in  the 
Lagrangian  description).  In  order  to  keep  this  paper  as 
simple  as  possible,  we  choose  the  Lagrangian  description 
and  concentrate  on  the  physical  nonlinearity,  which  seems  to 
be  the  most  significant  for  water  saturated  sandstones.  The 
assumption  of  a  periodic  microstructure  is  another  significant 
physical  assumption  made  for  the  purpose  of  simplifying  the 
effects  of  nonlinearity. 

We  consider  in  detail  a  layered  fluid/solid  medium.  The 
general  form  of  the  nonlinear  equations  is  derived,  but  we 
consider  only  linear  wave  solutions  in  this  paper.  Nonlinear 
waves  and  the  interplay  between  the  solid  and  fluid  nonlin¬ 
earity  will  be  explored  in  a  separate  paper.  The  stratified 
medium  is  perhaps  the  simplest  realization  of  a  porous  me¬ 
dium,  and  has  been  examined  in  several  studies,  both 
theoretical'^"*^  and  experimental.^  The  first  study  by  Rytov'^ 
considered  wave  motion  parallel  to  the  layers,  and  showed 
the  existence  of  two  waves:  the  fast  and  the  slow.  Bedford'^ 
later  showed  that  the  slow  wave  is  the  long-wavelength 
manifestation  of  the  second  mode  of  the  system.  This  is  a 
rich  system  for  wave  motion,  particularly  when  the  waves 
are  allowed  to  propagate  in  oblique  directions.  The  definitive 
theoretical  and  experimertal  studies  of  Schoenberg  and 
co-workers^’'^’'"^’*^  describe  the  wave  motion  and  slowness 
surfaces  for  fast  and  slow  waves.  We  give  for  the  first  time 
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TABLE  I.  Material  parameters,  from  Ref.  3. 


Material 

Compressional  speed 
(m/s) 

Shear  speed 
(m/s) 

Density 

(kg/m^) 

Aluminum 

6450 

3150 

2700 

Plexiglass 

2700 

1380 

1200 

Water 

1490 

1000 

the  complete  set  of  Biot  equations  for  the  stratified  fluid/solid 
system,  which  includes  fluid  viscosity  and  elastic  anisotropy. 

The  paper  proceeds  as  follows.  The  nonlinear  equations 
of  motion  are  outlined  and  the  scaling  procedure  is  defined  in 
Sec.  I.  The  main  difficulty  in  any  averaging  theory  of  hetero¬ 
geneous  media  is  to  calculate  the  “effective”  parameters  (see 
Table  I),  which  require  solving  simpler  problems  on  typical 
unit  cells.  Explicit  solutions  can  be  obtained  for  the  specific 
case  of  the  layered  medium,  and  are  addressed  in  Sec.  II.  All 
moduli  and  parameters  can  be  found  for  this  simple  geom¬ 
etry;  in  particular,  the  time-dependent  viscodynamic  operator 
is  explicit.  The  effective  nonlinear  equations  of  motion  are 
summarized  in  Sec.  III.  Section  IV  focuses  on  the  linear 
limit,  and  comparisons  are  made  with  Schoenberg’s 
model for  a  layered  medium  of  fluid  and  isotropic  solid 
in  Sec.  V. 

I.  GOVERNING  EQUATIONS  AND  SCALINGS 
A.  The  primitive  field  equations 

The  volume  V  comprises  fluid  and  solid  regions,  Vf  and 
,  with  boundary  dV  between  them.  The  displacement 
fields  are  w,-  and  1/^  in  the  solid  and  fluid,  and  the  stress  fields 
are  (Tj^  and  respectively.  All  equations  are  defined  in 
terms  of  the  reference  or  Lagrangian  coordinates.  The  gov¬ 
erning  equations  of  motion  and  constitutive  relations  are 


p,S^Ui=DjO-ji, 

(la) 

(lb) 

in  , 

Pfd^Ui  =  Djlji, 

(2a) 

Iji  =-p  Sji  +  2  v^,{Dj  Ui} , 

(2b) 

p=.AD,Ud, 

(2c) 

in  Vy,  and  the  continuity  conditions  on  dV  are 

(3a) 

nj1ji  =  nj(Tji. 

(3b) 

Here,  Fij  —  DjUf  is  the  deformation  gradient  tensor  in  the 
solid,  {aij}^(aij  +  aji)/2  -  Sija}^/^/3  is  the  symmetric  devia- 
toric  part  of  a  second-order  tensor,  and  rij  is  the  unit  normal 
on  SV  directed  into  the  fluid  region.  The  summation  of  re¬ 
peated  Latin  suffices  over  1,  2,  and  3  is  understood. 

Both  the  solid  and  fluid  constitutive  relations  are  nonlin¬ 
ear.  We  ignore,  however,  the  “geometrical”  nonlinearity  in 
the  fluid  equation  (2b)  for  the  sake  of  simplicity  (conceptu¬ 
ally,  many  of  our  conclusions  remain  valid  without  thi^  as¬ 
sumption).  Also,  the  equation  of  state  for  the  pressure,  Eq. 


(2c),  is  normally  expressed  in  terms  of  the  current  (Eulerian) 
density,  but  we  take  a  slightly  different  approach,  and  as¬ 
sume  that  it  depends  upon  the  dilatation  in  the  reference 
description.  This  is  simpler,  and  leads  to  the  same  nonlinear 
effects  that  are  normally  present  in  homogeneous  fluids  up  to 
second  order.  Also,  we  assume  that  the  fluid  equation  of  state 
possesses  an  inverse, 

(4) 

The  material  nonlinearity  of  the  phases  is  then  reflected  in 
the  nonlinear  behavior  of  the  functions  and  (jy,(F^„). 
The  stress  functions  aj^  and  the  pressure  are  normally  de¬ 
rived  from  strain  energy  potentials,  according  to 

dE,  _  dEf 

d{D;Ud' 

where  E^iDjU^)  and  Ejr(D^Ui)  are  the  elastic  strain  energy 
functions  for  the  solid  and  fluid,  respectively.  The  homogeni¬ 
zation  theory  outlined  below  is  not  dependent  upon  the  ex¬ 
istence  of  the  strain  energies,  and  the  final  equations  do  not 
involve  them.  However,  we  will  see  that  the  effective  po- 
roelastic  medium  also  possesses  a  strain  energy  function 
when  Es  and  Ej  exist. 

B.  Asymptotic  scaling 

In  order  to  simplify  the  system  of  equations  (l)-(3)  we 
employ  the  method  of  homogenization.  We  assume  the  me¬ 
dium  is  characterized  by  two  distinct  spatial  lengths, 
associated  with  the  micro-  and  macroscales,  respectively. 
The  number,  is  a  small  parameter,  6<^1.  Further¬ 

more,  we  assume  that  the  viscosity  scales  as 

r)=e^7j.  (6) 

References  7-9  and  1 1  provide  a  clear  motivation  for  this 
choice.  The  results  generated  by  the  homogenization  tech¬ 
nique  are  very  sensitive  to  the  presence  or  absence  of  the 
multiplier  ^  in  the  viscosity.  Without  it  the  “macro”  equa¬ 
tions  turn  out  to  be  of  viscoelastic  type,  whereas  the  presence 
of  the  term  ^  in  Eq.  (6)  leads  to  Biot-type  equations,  as  we 
will  demonstrate.  We  refer  the  reader  to  the  cited  papers  for 
further  details  of  the  homogenization  procedure.  [See  also 
the  comment  after  Eq.  (32).] 

We  assume  a  two-scale  expansion’^  in  terms  of  the  slow 
or  macroscopic  spatial  variable  x,  and  the  fast  variable 
y=x/6,  such  that  the  spatial  differential  operator  is 

Dj=d,+€-'dy..  (7) 

The  fast  variable  y  reflects  the  small-scale  structure  in  the 
problem  through  the  perturbation  parameter  e.  Our  goal  is  to 
eliminate  the  explicit  dependence  on  y,  leaving  us  with  equa¬ 
tions  in  X.  The  governing  equations  (l)-(3),  combined  with 
Eq.  (6),  become 

ePi + dy)  o-ji .  (8a) 

o-ji  =  (rji[(d.,^+  '  Syju„]  (8b) 

in  the  solid  phase,  and 

epfd^Ui=i€S,+dy)lji,  (9a) 
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• 

-P^ji+  f2  77r?,{(ed'^.+  dy.)Ui}, 

(9b) 

ap)  =  (S,+^~'dy)U^ 

(9c) 

in  the  fluid.  We  consider  the  ansatz 

M,(x,r,e)  =  u?(x,r)+S  e"u1{x,y,t), 

rt=  1 

ac 

/7(x,f,6)=p®(x,f)+ X 

n=  1 

X 

f/,(x,/,e)=2  e"t/?(x,y,r), 

n  =  0 

(10) 

X 

o-.j(x,r,e)  =  S  e"aU\,y,t), 

■'  n  =  0 

Substituting  these  expansions  into  the  governing  system  and 
comparing  like  powers  of  6  yields  a  sequence  of  asymptotic 
equations.  The  leading-order  equations,  of  order  unity,  are 

dy. 0-^  =  0, 

(11a) 

(lib) 

for  X  and  y  in  , 

ii 

o 

(12a) 

(12b) 

dy.l/!=0. 

(12c) 

for  X  and  y  in  Vf,  and  on  dV, 

(13a) 

nj^=njcr% 

The  next  set  of  equations,  of  order  e,  is 

a3b) 

Ps3^U^=^xjO-p+^yj^]i’ 

(14a) 

,  ..2^ 

^  Jl  nr  ^  <i  yp 

or  pq 

(14b) 

in  Vj, 

Pfd]U^,=d,%%+dyf.],, 

(15a) 

Vjr-p'Sji+lvdA^V'l)^ 

asb) 

0{p^)  =  d,V^,+dyV\, 

(15c) 

in  Vf,  and 

U]  =  u\, 

(16a) 

nj^ji^njCTji, 

(16b) 

on  dV. 


C.  Homogenization  of  the  asymptotic  equations 

Introduce  the  relative  fluid  displacement  vector 
defined  as 

w®(x,y,r)  =  <^[l/°(x,y,t)-M”(x,r)],  (17) 

where  0<<f><  1  is  the  porosity  or  volume  fraction  of  the  fluid, 
VjfV.  We  can  then  rewrite  Eq.  (15c)  as 

The  dependence  upon  the  fast  scale  may  now  be  eliminated 
by  averaging  over  y.  A  formal  definition  of  the  averaging 
procedure  can  be  given  for  nonperiodic,  statistically  defined 
media;  see,  for  example,  Burridge  and  Keller.  Here,  for  the 
sake  of  simplicity,  we  assume  periodicity  on  the  small  scale, 
so  that  the  average  is  trivial.  Thus  integrating  Eq.  (18)  over 
Vf,  using  the  displacement  continuity  conditions  (16a)  and 
the  assumed  periodicity  in  y,  gives 

4>{  - ap°) + + yj ' 

=  (l-(^)V7‘f  ^iV(y),  X  in  V;:.  (19) 

^  s 

Substitution  of  Eqs.  (12b)  and  (15b)  into  Eq.  (15a)  yields 
Pfd^w°  +  dy,{p  ‘  (f>Sji-2  rj^A^yy^Vs) 

=  -<i>[Pfdyi^d,p%  x,y  in  Vf,  (20) 

while  Eqs.  (12c)  and  (13a)  become,  respectively, 

(?  w°=0,  y  in  Vf,  (21a) 

w®=0,  y  on  dV.  (21b) 

Next,  integrating  Eq.  (14a)  over  V^  and  Eq.  (15a)  over  V^, 
then  adding  the  results  and  using  the  continuity  conditions 
(16b)  and  periodicity,  we  find 

pd]u%\,i)  +  PfV]'  [  d]w^i  dV{y) 

JVf 

= -  -  ^)y7'  ^ 


where 

p=(l>Pf+{l-(l>)Ps 

is  the  average  density.  Finally,  using  Eq.  (12b),  we  rewrite 


the  boundary  condition  (13b)  as 


rijCrl^-p^rii 


(24) 


The  system  of  equations  (20)  and  (21)  forms  a  closed 
boundary-value  problem  with  respect  to  ^;W^(x,y,0»  ^nd  in 
principle,  one  can  express  this  function  in  terms  of  w,-  and 
p^.  Similarly,  Eqs.  (11)  and  (13b)  give  a  well  defined 
boundary- value  problem  with  respect  to  W  / ,  which  envies  us 
to  present  m-  as  a  function  or  functional  of  and  /?  .  Both 
problems  are  solved  explicitly  for  the  case  of  a  periodically 
layered  medium  in  the  next  section.  Inserting  the  results  in 
Eqs.  (19)  and  (22),  we  arrive  at  a  closed  master  system  of 
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equations  with  respect  to  and  w^  ,  with  the  averaged 

relative  fluid  displacement  defined  as 

w,(x,f)  =  V7'  [  ^v®(x,y,^)^/V(y).  (25) 

}Vf 

II.  CONSTITUTIVE  THEORY  FOR  A  LAYERED  MEDIUM 

We  consider  periodically  alternating  fluid  and  solid  lay¬ 
ers.  The  thickness  of  each  fluid  layer  is  /  in  terms  of  x  and 
L=l/e  in  terms  of  y.  Let  3;=y*ii  be  the  coordinate  in  the 
direction  of  layering.  The  two  microproblems  defined  at  the 
end  of  the  previous  section  then  depend  only  upon  y,  and 
hence  reduce  to  unidimensional  problems.  This  is  the  great 
simplification  that  arises  from  the  layering,  and  it  is  not 
present  for  any  other  configuration. 

A.  The  permeability  operator 

The  permeability  operator  results  from  the  solution  of 
Eqs.  (20)  and  (21).  Let 

Vi=d,w^,  =  (26) 

then  Eqs.  (20)  and  (21b)  become 

dfV dyp  ‘  +  7)pJ '  d]vi  +  <^//( X, t ) ,  (27a) 

nidyVi  =  Q.  (27b) 

We  consider  y -periodic  solutions  of  this  system.  The  bound¬ 
ary  conditions  are  those  for  a  rigid  wall  at  y  =  0  and  y  =  L, 
which  combined  with  the  incompressibility  condition  for  wf , 
implies 

niVi  =  ^<^v=PijVj,  (28) 

where  P(n)=I-n(S)n  projects  on  to  the  horizontal  plane. 
Contracting  Eq.  (27a)  with  ri  and  using  Eq.  (27b)  allows  us 
to  eliminate  the  pressure  gradient  term.  Bearing  in  mind  Eq. 
(28)  we  deduce  that  Vi  satisfies 

d^Vj- vd]vi=(}>Pijfj{^J).  (29) 

where  v=  pip j  is  the  kinernatic  viscosity.  This  can  be  solved 
by  standard  means,  i.e.,  using  a  Fourier  series  in  y.  A  very 
similar  type  of  problem  is  discussed  by  Sneddon.  Taking 

into  account  the  initial  data  (i;  =  0  for  r<0),  we  find 

^2  nirv 

Vi(x,ya)  =  (PPijZ  -  [  1  -  ( -  1 )”] sin 

n=l  n  L 

X  I'  di  (30) 

Jo 

Averaging  y,(x,y,r)  over  the  layer  yields  i;/(x,r),  which  sat¬ 
isfies 

Kij{t)  =  PijK{t),  (31) 

where  *  denotes  convolution,  and 

'TT  n  =  1  ( i  /  \  ^  1 


FIG.  1.  The  permeability  function  for  the  stratified  medium  in  nondimen- 
sional  units. 


Hence  the  permeability  of  the  medium  is  a  projection  opera¬ 
tor  onto  the  transverse  plane.  Also,  K{t)  is  independent  of  e, 
because  vllJ—pIpjP'. 

Note  that  Eq.  (31)  implies 

=  +  (33) 

where  K{0)  =  (j)  fox  the  layered  medium,  and  K\t)  is  plotted 
in  Fig.  1.  It  follows  from  Eq.  (32)  as 

Sv  ^  I  ^  ^ 

K'{t)=-(f>jj  expl -7r^(2n- 1)^  pj.  (34) 


This  is  a  well  behaved  and  convergent  series,  except  for 
values  of  t  approaching  zero,  where  it  is  evidently  singular. 
The  behavior  near  zero  can  be  seen  by  transforming  the  sum 
using  the  Poisson  summation  formula:  for  a  given  function 
g{x). 


00 

2  gin) 

n=~x 


-  il  TTtnx 


dx. 


(35) 


Applying  this  to  the  sum  in  Eq.  (34)  and  performing  the 
integration  yields 


K'{t)=-c}> 


-+2  (-1)"  exp 

2  «=l 


n^L^\ 
4vt  I 


(36) 


The  function  K'{t)  therefore  has  an  integrable  t  singu¬ 
larity  at  zero. 


B.  The  stress  and  strain  in  the  solid  layers 

We  now  turn  to  the  system  defining  the  displacement  w  ■ 
in  terms  of  the  macroscopic  parameters  w?  and  p®,  Eqs.  (11) 
and  (24).  Proceeding  as  we  did  for  the  previous  microprob¬ 
lem,  we  look  for  the  affine  solution  depending  on  the  vertical 
coordinate  y  only.  Hence 

(37) 

where  the  vector  does  not  depend  upon  y.  Then  u]  auto¬ 
matically  satisfies  the  equilibrium  equations  within  the  solid 
region,  Eq.  (11).  The  traction  boundary  conditions  of  Eq. 
(24)  become,  using  Eq.  (11b), 
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nj(rji{d^^u°  +  n„G„)=-p°ni.  (38) 

This  is  an  algebraic  system  of  three  equations  in  three  un- 
knowns,  with  solution 


The  effective  strain  energy  follows  from  Eq.  (Al),  which 
becomes  in  the  present  notation 

1V(F,0  =  ( 1  -  <^)£,[F.,  +  n,G/F,p)]+  ^Ej{0{p)). 

(45) 


G„=G„(,?,«V). 


(39) 


Explicit  versions  of  this  will  be  discussed  later,  but  for  the 
moment  the  implicit  solution  of  Eq.  (39)  is  sufficient.  The 
leading-order  stress  in  the  solid  then  follows  from  Eqs.  (11b) 
and  (37)  as 


O’/I  =  + MmG„ ) = . 


(40) 


III.  MACROSCOPIC  EQUATIONS 
A.  The  general  nonlinear  equations 

We  are  now  in  a  position  to  state  the  leading-order  mac¬ 
roscopic  equations  of  motion.  The  fundamental  field  vari¬ 
ables  are  the  solid  displacement  M?(x,r)  and  the  averaged 
fluid  displacement  w,(x,r),  and  their  governing  dynamic 
equations  follow  from  Eqs.  (22),  (26),  and  (31).  In  order  to 
simplify  the  notation,  we  make  the  replacements 
M?(x,r)-^M,(x,f)  and  p°(x,f)— p(x,/).  The  dynamic  equations 
are 


p^,Ui{\,t)+Pfd]wi(,x,t)=Tjij{\,t),  (41a) 

d,Wi{x,t)=-Kij(t)*{d]uj(x,t)  +  p]'p^j{\,t)'\,  (41b) 

where  is  given  explicitly  in  Eq.  (32).  The  stress  and 

pressure  are  related  to  the  displacements  by  the  right-hand 
side  of  Eq.  (22)  and  Eq.  (40),  and  by  Eqs.  (19)  and  (37), 
respectively,  or 

Tj-,(x,r)  =  ( 1  -  <^)rj,(F,p)-  .^.p(x,t)  <5,, ,  (42a) 

^(x,r)  =  <Ae(x,r)  -  <^<?(p)-(l  “</»)«, ■G,(F,p),  (42b) 

where  =  d^  ui  =  m/j  is  the  macroscopic  deformation 
gradient  tensor,  e  =  Ui  iy  and  —  is  the  relative  fluid 
dilatation.^ 

Equations  (42)  provide  the  bulk  stress  Tj^  and  the  rela¬ 
tive  fluid  dilatation  ^  in  terms  of  the  solid  deformations  F ij 
and  the  fluid  pressure  p.  When  the  pressure  is  zero,  the  bulk 
stress  is 


t^,.  =  (1-<A)7’;,(F,0),  p^O.  (43) 

Hence  we  can  identify  (1  “<^)7'y,(F,0)  as  the  stress  function 
for  the  dry  frame,  or  the  “open”  system.^  If  the  pores  are 
sealed,  or  “closed,”^  then  the  associated  stress  function  fol¬ 
lows  from  Eqs.  (42)  with  ^=0. 


B.  Energy  potentials 

The  macroscopic  elastic  constitutive  relation  for  the 
stress  in  Eq.  (42a)  can  be  related  to  the  effective  strain  en¬ 
ergy.  Comparison  of  Eqs.  (42)  and  (A  13)  implies  that  the 
former  can  be  replaced  by  the  simpler  relations 


dW  _  dW 


(44) 


Hence  the  effective  medium  is  elastic  in  the  sense  that  it 
possesses  a  stored  energy  potential  W.  Note  that  W  is  a  func¬ 
tion  of  the  kinematic  strains  F  and  but  its  functional  form 
is  defined  by  F  and  p.  The  remaining  identity  (42b)  is  an 
implicit  relation  for  p  in  terms  of  F  and  i.e.,  p=p(F,f). 

The  related  potential  11  is  defined  in  the  Appendix.  It 
follows  from  Eq.  (A7)  that  the  two  constitutive  relations  of 
Eq.  (42)  are  equivalent  to 


<?n 


(46) 


The  functional  definition  of  n=n(F,p)  is  given  in  Eq.  (A6), 
which  becomes  in  the  present  notation 


n(F,p)  =  (l-<^)£,(F^,  +  «,G/F.p)) 


-b  ( 1  -  <^)p«mG„(F,p)  -  4>pe 


+  <)>[  av)dv-  (47) 

This  is  an  explicit  equation  for  the  potential,  apart  from  the 
fact  that  the  functions  G,(F,p)  are  given  by  the  implicit  re¬ 
lation  (38).  A  simpler  form  for  11  can  be  obtained  by  noting 
that  it  is  a  partial  Legendre  transform  of  the  energy  function 
1V(F,^).  The  differential  relations  (44)  and  (46)  imply  the 
connection 

w=n+^p.  (48) 


IV.  LINEAR  THEORY  FOR  THE  STRATIFIED  MEDIUM 

It  is  instructive  to  consider  the  limiting  case  for  small 
amplitude  waves,  for  which  the  linear  constitutive  equations 
in  the  solid  and  fluid  are  sufficient.  The  individual  energy 
potentials  are  now 

=  Ef{e)  =  \Kfe‘^,  (49) 

where  Kf  is  the  fluid  bulk  modulus  and  the  elastic  moduli 
possess  the  symmetries  Cij^=Ci^iij  and  The 

linear  stress  in  the  solid  is 

o'jiiu„,„)  =  CjiuUk.i-  (^®^ 

The  deformation  gradient  tensor  F  can  be  replaced  in  all 
expressions  by  the  macroscopic  linear  strain  tensor 

e=(F+F^)/2,  or  Equation  (38)  becomes 

2  Mnear  system  for  G/  which  can  be  solved  easily,  yielding 

G,(e,p)=-fc,*,e*,-e,V«tP,  (51) 

where 

Qik-^ijkl^ ^ikl^Qip  ^pjkl^j  • 

The  functions  Tji  of  Eq.  (40)  then  follow  as 
^ ;/ ( j^nklQmn ) ^ /c/ ^kij^ kP * 


(52) 


(53) 


Using  the  fact  that  (2{p)—  —piKf,  whiwh  follows  from  Eqs. 
(5)  and  (49),  the  constitutive  relations  of  Eq.  (42)  become 
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Tji{xj)  =  Cl„e„-M-'Mijp,  (54a) 

a^,t)  =  M-'Mijeij  +  M-'p,  (54b) 

where 

( 1  -  (f>)iCijkrbmijb„kiQmn),  (55a) 

M  =  [(f)/ Kf+{\  — 4>)ninjQij^]  (55b) 

Mij  =  [<t>8ij+  ( 1  -  (55c) 


Note  that  C°  are  the  dry  frame  or  open  pore  moduli; 
see  Eq.  (43).  Equations  (54)  are  equivalent  to 


(56a) 

p  =  MC-Mijeij,  (56b) 

where 

^‘ijkl~  ^ijkl  ^  ^ij^kl '  (5"^) 

are  the  “closed  pore”  moduli  (^=0).  Also,  (l-<^)  'C®*/ 
can  be  identified  as  the  moduli  for  a  state  of  plane  stress  in 
the  elastic  solid.  Thus  Eq.  (54a)  indicates  that  macroscopic 
plane  stress  prevails  when  the  pressure  is  zero,  as  expected. 

The  potential  11  follows  by  integrating  r,,  and  C,  using 
Eqs.  (46).  The  energy  density  then  follows  from  Eq.  (48).  We 
find,  after  some  simplification,  that 


n  2  ^ ijkl^ij^kl 


The  former  is  a  remarkably  simple  and  physically  appealing 
form.  It  clearly  separates  the  strain  energy  into  solid  and 
fluid  parts,  where  the  latter  depends  upon  the  effective  bulk 
modulus  M. 


V.  ISOTROPIC  ELASTIC  LAYERS  AND  WAVES 
A.  The  general  equations 

Let  \  and  p.  be  the  Lame  moduli  of  the  solid;  then  Eqs. 
(52)  become 

Q-'  =  (X  +  2/t)-'n®n+/i.-‘P(n),  (59a) 

X 


biki= 


X.  +  2/li 


(59b) 


It  is  a  simple  matter  of  algebra  to  show  that  the  parameters  in 
Eqs.  (55)  are  now 

(60a) 
(60b) 


ijP  ki'^2plijkl), 

X  +  2/1/ 


M  =  l  —  + 

.Kf 


Mir 


riiTij- 


X  +  2/x<^ 
X  +  2yLt 


Pij 


M, 


(60c) 


where  Ijki^iP ikP ji  + P uP ikV'^  is  the  “in-plane”  fourth- 
order  identity  tensor  and  Xo=2/iX/(X+2/i.)  is  the  plane 
stress”  Lame  modulus.  The  two  constitutive  equations  are 
(54a)  and  (54b),  the  former  simplifying  to 


r„= ( 1  -  <^)(XoeP,v+ 2/te,;)  -M" 


(61) 


where  erhiki^ki  in-plane  strain,  and  €  =  ejj  .  The 

constitutive  relations  (54b)  and  (61)  are  written  with  the 
solid  strain  and  the  fluid  pressure  p  as  the  fundamental 
variables  for  the  poroelastic  medium.  Alternatively,  if  we 
choose  the  strains  efj  and  ^  as  the  primitive  variables,  then 
the  constitutive  relations  are  Eqs.  (56).  The  energy  potential 
W  of  Eq.  (58)  reduces  to 

1 

w=-{l-<f>){  Xo^^  +  Ipiijeij)  +  — .  (62) 

The  Biot  equations  are  those  of  a  transversely  isotropic  po- 
roelastic  medium,  as  expected.  The  general  Biot  equations 
for  a  material  with  this  symmetry^®  have  eight  independent 
moduli:  five  for  the  solid,  one  for  the  fluid,  and  two  for  the 
fluid/solid  interaction.  It  is  interesting  that  the  present  equa¬ 
tions  only  have  five  independent  moduli,  two  each  for 
and  Mij,  and  one  for  M.  Three  of  the  eight  possible  con¬ 
stants  are  identically  zero,  which  can  be  ascribed  to  the 
plane-stress  configuration  of  the  solid. 

The  equations  of  motion  can  be  further  simplified.  Sub¬ 
stituting  Tij  from  Eq.  (61)  into  Eq.  (41a)  gives  a  decoupled 
system  of  equations.  We  now  let  n=e3,  where  (61,62,63) 
form  an  orthonormal  triad.  Then,  using  Eq.  (33),  we  find 

p<?^W3+p  3  =  0, 

=  (l“<^)"*A^'(r)*(P/^?Wa  +  P.a)» 

where  the  Greek  suffices  a  and  /S  are  restricted  to  I  and  2, 
and  the  summation  convention  is  implicit.  Finally,  eliminat¬ 
ing  C  between  Eqs.  (41b)  and  (56b)  gives 

where  K  is  given  in  Eq.  (32).  The  four  equations  in  Eqs.  (63) 
and  (64)  form  a  closed  set  of  equations  of  motion  for  u(x,r) 
and  p(x,r).  Next,  we  consider  four  examples  of  wave  solu¬ 
tions  in  the  layered  medium. 

B.  Examples  of  wave  motion 

1.  Propagation  in  the  vertical  direction 

As  a  first  example  consider  motion  in  the  vertical  direc¬ 
tion  only:  u=M(jC3,Oe3  and  w=w(x3,f)e3.  Then  Eq.  (41b) 
implies  w  =  0  while  Eq.  (42)  gives  Tjre23Mji  and 
The  equation  of  motion  (63a)  becomes  simply 

pd^u(xj)-Mdlu{xj)^0.  (65) 

This  is  a  scalar  wave  equation  with  wave  speed  ^jM/p, 

2.  A  damped  shear  wave 

Consider  shear  motion  in  the  horizontal  plane: 
,r)e2  and  w=w(xi,r)e2.  The  dynamic  equations 
again  reduce  to  a  single  equation  for  u : 

P,d^u{x,t)- =  <pr'K'it)*Pfd^u{x,t). 

(66) 

The  shear  wave  propagates  with  a  time  delayed  memory 
function  for  the  inertial  term.  The  appearance  of  K'  reflects 
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the  influence  of  the  fluid  viscosity  on  the  shear  wave  in  the 
solid. 


3.  Wave  motion  when  the  fluid  is  inviscid: 
Schoenberg's  problem 

The  limiting  case  of  an  inviscid  fluid  is  relevant  to  ex¬ 
periments  performed  by  Fiona  et  al?  The  fast  and  slow 
waves  in  the  system  are  both  nondispersive  and  nonattenuat¬ 
ing  in  this  limit.  But  the  anisotropy  of  the  configuration 
means  that  both  wave  types  exhibit  directional  dependence, 
which  can  be  best  understood  by  considering  the  slowness 
surfaces.  Schoenberg derived  the  equations  for  the  slow¬ 
ness  surface  from  the  dispersion  relation  for  the  long- 
wavelength  modes  of  a  layered  system  of  fluid  and  isotropic 
solid.  We  will  demonstrate  that  exactly  the  same  slowness 
surface  follows  from  the  governing  Biot  equations  outlined 
above. 

The  inviscid  nature  of  the  fluid  means  that  the  operator 
K  is  instantaneous,  or  /l'  =  0.  Differentiating  Eq.  (64)  with 
respect  to  time  then  implies,  using  Eq.  (33),  that 

+  MijCij)  =  MK{Q){d]i+p]  'p,aa).  (67) 

where  A'(0)  =  (/>.  Consider  motion  in  the  ^^1-^:3  plane  with 
solutions  of  the  form 


Wi(/-5iXj  -53X3),  W9  =  0,  w3(r*-5iXi  -53^3), 

P=p{t-SyX^-S^X2,). 

where  and  53  are  the  horizontal  and  vertical  components 
of  the  slowness.  Let  v=^^u;  then  Eqs,  (63)  imply 

pUs”  53P  =  0,  (69a) 


( 1  -  cl,s])v ,  -  Vl  -  clilclstp  =  0,  (69b) 

where  Cp  and  Cpi  are  the  wave  speeds  for  longitudinal  waves 
in  the  solid  in  bulk  and  as  a  thin  plate,  respectively.  That  is. 


c;={k  +  2fi)/p,,  c^,  =  (Xo  +  2m)/p,,  (70) 

and  therefore  Ko/2p.  =  \/l  Eliminating  V'i  be- 

tween  Eqs.  (67)  and  (69a)  yields  the  third  equation, 


(l“^)Vl-c“//c;5,i;i- 


p  =  0,  (71) 


Equations  (69b)  and  (71)  are  satisfied  if  Si  and  53  are 
related  according  to 


P 


^  {  -2  2\j_ 

-(Cf  -5,)  + 


(l-«^) 

p,{\-clis]) 


where  Cf  =  \jKf/pf  is  the  speed  of  sound  in  the  fluid.  Equa¬ 
tion  (72)  defines  the  slowness  surface  for  the  fast  and  slow 
waves,  and  it  agrees  with  Schoenberg’s*^  equation,  derived 
by  taking  the  low-frequency  limit  of  the  exact  dispersion 
relation  for  a  finely  laminated  medium  of  alternating  solid 
and  fluid  layers  (see  also  Refs.  3  and  14).  We  may  rewrite 
Eq.  (72)  as 


0p 


2\  (“^fast  •^l)(*^slow  *^1)’ 


(73) 


FIG.  2.  The  slowness  surfaces  for  two  different  stratified  media,  consisting 
of  water/aluminum  (Al)  and  water/plexiglas  (PI).  The  material  parameters 
are  in  Table  I  and  the  porosity  is  (^=0.5  in  each  case.  The  ordinate  and 
abscissa  are  the  dimensionless  horizontal  and  vertical  slownesses  s, 
=  j|  \j<t>M/pf  and  Sj  =  respectively.  The  ellipse  of  Eq.  (76), 

which  becomes  a  circle  in  these  units,  is  also  shown.  The  plate  slowness  c"/ 
corresponds  to  s,  =0.268  for  aluminum  and  s,  =  0.561  for  plexiglas. 


where  ifast<'Ssiow  3^^  horizontal  slowness  for  the  fast  and 
slow  waves  propagating  horizontally.  They  satisfy 


Pf 


('^fast'*’'^slow)" 


<!> 


<f>  l-4>\  1 

- !■ - 

Pf  Ps  I  Cp, 


(74a) 


^22  ^ 
^  ‘^fasrslow 

Pf 


<f>  1 

- _+ - ^  _ 

Pf^f  Ps^p  I  ^ pi 


(74b) 


These  slowness  values  were  also  obtained  by  Rytov*^  who 
analyzed  waves  traveling  horizontally  in  a  fluid/solid  layered 
medium. 

Schoenberg*^  considered  a  Plexiglas/water  system  with 
<^=0.6,  while  Fiona  et  al?  considered  the  same  system  with 
4>=0.5,  and  also  aluminum/water  with  <^=0.5.  In  all  cases  it 
was  found  that  therefore  the  slowness 

surface  of  Eq.  (73)  comprises  an  inner  closed  sheet  for 
corresponding  to  the  fast  mode,  and  an  outer 
sheet  for  which  gives  the  slow  mode.  A  stop 

band  exists  for  •  The  value  of  is  very  close 

to  Cpi  in  all  these  cases,  and  very  accurate  approximations 
can  be  obtained  by  iterating  from  this  starting  value;  thus 


2  1 

■^slow~  2  + 

<f>PsCl 

_  Pf 

(75a) 

2  ^  _ 

(75b) 

■*fast  2 

^pl 

<f>Pscll 

Thus  the  slowness  surface  is  essentially  the  ellipse 

I  $2,  1 

Pf  p  M' 


(76) 


punctuated  by  a  small  stop  band  to  the  left  of  =  Cp/*.  The 
slowness  surfaces  for  two  material  combinations  are  plotted 
in  Fig.  2,  which  also  shows  the  ellipse  of  Eq.  (76). 

The  displacement  polarization  follows  from  Eqs.  (69a). 
For  a  given  value  of  ,  we  have 


U|  _  S,  yjl-cjtlcl 
P  Ps{l-cl,s])'  p  p’ 


(77) 
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where  53  is  defined  by  Eq.  (73).  Assume  that  5i>0  and 
53>0,  as  in  Fig.  2.  The  ratios  v^lp  and  v^ip  are  both  posi¬ 
tive  for  the  fast  wave  (the  branch  The  same  ratios 

are  negative  and  positive,  respectively,  for  the  slow  wave. 
Thus  the  horizontal  velocity  is  in  phase  with  the  acoustic 
fluid  pressure  for  the  fast  wave,  but  out  of  phase  for  the  slow 
wave. 


4.  Dynamic  compatibility 

The  phenomenon  of  dynamic  compatibility  was  noted 
by  Biot  in  his  first  paper  on  the  dynamics  of  isotropic  porous 
media.  ^  If  the  material  parameters  satisfy  a  certain  constraint, 
then  a  wave  solution  exists  which  has  no  viscous  attenuation 
and  its  relative  fluid  motion  is  zero,  or  w=0.  Analogous 
compatibility  conditions  exist  for  the  layered  medium,  as  we 
now  demonstrate  for  horizontal  wave  motion.  Consider  lon¬ 
gitudinal  motion  in  the  horizontal  plane,  u=w(;ri,r)ei  and 
w=w{x^  ,t)ej.  This  leads  to  a  pair  of  coupled  equations  for 
u{x,t)  and  >v(x,r). 


pd^u-[(\-<t>){\()'^2/x)  +  Mq^]dlu  +  pfd^w 
—  Mqd^w  —  O, 

djW  +  K{t)* 


dju - dZ{qu~\-w) 

Pf 


=  0, 


(78a) 

(78b) 


where  ^  =  (X  +  2/i,<;6)/(X  +  2/x).  The  system  reduces  to 
Pfd]u-Mqdlu-=^.  w  =  0,  (79) 


if  the  physical  parameters  are  related  by 
Mqlpf={{  1  <^)(Xo  +  2/i)  +  Af^“]/p, 

or  equivalently,  if  the  densities  are  in  the  ratio 
Xq  2yu<  2fx 

Xo(  1 ^)  +  2/Lt<^  M 


(80) 


(81) 


This  is  the  condition  of  dynamic  compatibility,  and  when  it  is 
satisfied  waves  can  propagate  unattenuated  by  viscous  drag 
with  wave  speed  yjMq/pf.  Alternatively,  the  right-hand  side 
of  Eq.  (81)  may  be  written  in  terms  of  the  bulk  compres- 
sional  sound  speeds  and  Poisson’s  ratio  of  the  solid, 
v^  =  k/2{\  +  fi),  yielding 


p.  - 

The  right-hand  member  is  an  increasing  function  of  (f>  for 
most  material  combinations,  i.e.,  Cp/c^(l  —  v^)- 1  >  1  nor¬ 
mally.  Therefore,  for  a  given  pair  of  materials,  there  is  a 
unique  value  of  the  porosity  at  which  compatibility  is  at¬ 
tained  as  long  as  pjpf>{\-  vj/  . 


VI.  CONCLUSION 

The  theory  of  homogenization  has  been  used  to  derive 
the  nonlinear  equations  for  a  fluid  infiltrated  solid  skeleton 
starting  from  the  fundamental  equations  governing  the  mo¬ 
tion  on  the  microscale  of  the  pores  and  grains.  The  canonical 
nicrogeometry  of  a  layered  system  of  alternating  solid  and 
fluid  constituents  has  been  analyzed  in  detail.  The  main  re¬ 


sults  are  summarized  in  the  “averaged”  or  macroscopic 
equations  of  motion  (41),  and  the  nonlinear  constitutive  re¬ 
lations  of  Eqs.  (42).  The  viscous  effects  of  the  pore  fluid  are 
contained  in  the  permeability  convolution  operator  defined 
by  Kij{t)  of  Eq.  (32).  The  stress-strain  relations  in  Eqs.  (42) 
allow  for  arbitrary  nonlinear  behavior  in  the  elastic  solid  and 
in  the  fluid.  The  key  quantity  is  the  vector  function  G,  de¬ 
fined  by  Eqs.  (38)  and  (39).  We  have  also  shown  how  the 
stress  and  strain  are  related  to  a  macroscopic  potential  en¬ 
ergy,  in  the  spirit  of  Biot’s  later  work  on  nonlinear  mechanics 
of  porous  media.^' 

The  linear  limit  of  these  equations  is  of  interest,  and  we 
have  derived  for  the  first  time  in  Secs.  IV  and  V  the  full  set 
of  linear  Biot  equations  for  a  layered  medium.  The  predicted 
form  of  the  slowness  surface  for  isotropic  elastic  layers 
agrees  with  that  obtained  by  Schoenberg'^’*^  from  the  exact 
dispersion  relation  for  the  system.  This  gives  us  confidence 
in  the  general  validity  of  the  nonlinear  Biot  theory  as  derived 
here. 
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APPENDIX:  ENERGY  CONSIDERATIONS 

The  leading-order  term  for  the  macroscopic  energy  den¬ 
sity,  W,  in  the  porous  medium  is  a  linear  combination  of  the 
energy  in  the  solid  and  the  energy  in  the  fluid,  each  averaged 
over  its  domain.  The  fluid  energy  follows  from  Eq.  (15c)  as 
E^O{p^)\  which  is  independent  of  y,  so  that  the  sum  of  the 
average  partial  energy  is 

W=V~'  f  E,{d^  ul+d  u\)dV(y)  +  (l)Ef[0{p^)\ 

j 

(Al) 

The  purpose  of  this  appendix  is  to  convince  the  reader  that  W 
of  Eq.  (Al),  considered  as  a  function  of  the  independent 
macroscopic  strain  variables,  does  indeed  yield  the  stress  and 
pressure  as  partial  derivatives  according  to  Eq.  (44).  We  fol¬ 
low  Sanchez-Palencia^  and  first  consider  the  function 

V-'  I  EAd.jil+dyji^dV^y) 

+pV-' f  ujnidSiy).  (A2) 

The  point  is  that  A  is  clearly  a  function  of  the  arguments 
shown,  because  the  field  variables  w-  are  implicit  functions 
of  both  dr  and  p^,  as  discussed  above  [see  Eqs.  (11)  and 

m  “ 

(13b)].  The  first  variation  is 

Jv,  <?r,; 

X  {d,.Suyd,.Su'j)dV(y) 
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V '  f  ( d^.Su]  +  dy.Su] ) <T°j  dV( y) 


+  Sp°V-'\  ulnidS(y)  +  p°V-'\  Sujn^dSiy). 

(A3) 

The  terms  involving  Suj  annihilate  one  another  on  account 
of  Eqs.  (12b),  (13b),  and  the  first  of  (5).  We  are  then  left  with 


* 

M  =  dV(y) 


+  Sp'^V 


0,7-1 


/v 


^V(y). 


(A4) 


Thus,  referring  to  Eqs.  (37)  and  (40),  we  have 


n  7,  o>A 

( 1  -  «/.)n,G,(^^^«”,p‘')  = 

dA 


respectively. 

Consider  the  related  function 


uid,  uip<^)=A(d,y„,p^)-(f>p'^d,M° 


+  </>P 

Pzinh 


(A5) 


(A6) 


where  p^j^b  ambient  pressure.  H,  like  A,  is  a  function 
of  the  independent  variables  and  Taking  its  partial 
derivatives  with  respect  to  these,  and  using  the  identities  of 
Eqs.  (A5)  and  (42),  we  obtain  the  identities  of  Eq.  (46).  Thus 


Sn{V,p)  =  rjiSFj-^Sp 

(A7) 

or 

SmF,C)^TjtSFji+pS^, 

(A8) 

where 

ft(F,^)=n(F,p)+^/7. 

(A9) 

We  will  now  demonstrate  that  ft  and  W  are  one  and  the 
same.  The  latter  can  be  rewritten  using  Eq.  (15c)  as 


W=  V- '  E,{d,y+dyy)dV{y) 

+  Ef{d,Ul  +  dyUl)dV{,y).  (AlO) 

The  first  variation  is 

<5W=  V" '  [  ( d^. Su]  +  dy^j ) (T%  dV{y) 

-V-’f  p°S(^p°)dV{y),  (All) 

JVf 

where  Eqs.  (15c)  and  the  second  of  (5)  have  been  used.  The 
integral  over  the  fluid  region  can  be  simplified  further  using 
Eqs.  (16a),  (18),  and  integration  by  parts,  giving 


-<f>p^s{dy)-p<^s{d,w!) 

+p°V'  f  SulrtidSiy).  (A12) 

JdV, 

The  terms  with  Suj  again  disappear  because  of  Eqs.  (12b), 
(13b),  and  the  first  of  (5),  leaving 

5W=<!>(5,^M^|  V-'P  <rl„  dV{y)-(f>p^S„„^ 

-p°Sid,.Wi)-  (A13) 

This  implies  the  identities  of  Eq.  (44),  and  the  equivalence 
ft=W.  (A14) 
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